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ABSTRACT 


A computer program that can determine the collapse load 
of a transversely ioaded grillage that allows for various 
types of loading and failure modes 1s presented. The 
program allows the user to constrain particular grillage 
members to have fully developed plastic hinges at the 
location of maximum moment which gives the program the 
capability to study redundancy of transversely loaded 
Grillages. The lower bound theorem of Limit analysis and 
linear programming techniques are used to determine a grid 
load factor and provide a failure mode and associated 
bending moments. A user’s guide and program explanation are 
orovided. An example grillage is studied using this pregram 
and the load factor and redundancy is discussed. 
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INTRODUCTION 


The capabilility to easily and effectively study a 
structure from the stand point of redundancy would 
Significantly contribute to the tools available to engineers 
for calculating safety factors and provide requiatory bodies 
an inexpensive way to verify or develop safety factors with 
analytical basis. 

Strengthening of ship structures has been an area of a 
Significant amount ot work. Ice loading is particularly 
threatening to ship safety and also has the characteristic 
of very large localized or point loads. The program 
presented here is ideal for the study of local loading or 
failure. 

Previous work done on grillages that used similar 
techniques to determine the load carrying capability and 
possible failure modes was done by Tait and Hodge (1) in 
their develapment of a computer program called BEAMPLT that 
utilizes linear pragraming techniques to find a solution to 
transversely loaded grillages. [This work was later expanded 
upon to provide a simple method of computing failure modes 
and load factors under various patch loads by Abbott (2). 

Hadge (1). in his paper on transversely loaded 
Grillages. demonstrated the feasibility of using linear 


Programing to calculate a failure mode and load carrying 








capacity (load factor) for transversely loaded grillages. 
His work allowed for several types of boundary conditions. 
These conditions included simply supported. clamped. and 
cantilevered. The loading scheme was limited to point loads 
on each node. 

Abbott (2) used this approach to solve for various 
types of patch loads on a grillage. He wrote a computer 
program that employed the same techniques as Hodge used. 

He employed linear programing techniques for calculating a 
failure mode and load factor for a given grillage. The 
patch loads that he was interested in were those that could 
be associated with relatively small area ice loading on the 
side of a ship. These snatch loads were transferred to the 
nodes of the qrillage using a simple lever principle. 

Reth Hodge (1) and Abbott (2) allowed failures to occur 
only at nades in the grillage. Hodge’s scheme allowed 
laading anly at the nodes. Abbott’s method was adequate to 
estimate failure with reasonable accuracy, But his 
assumptions in transferring patch loads into nodal loads was 
simplistic and did not allow for point loads off of a node. 
Abbott*s approach followed Hodges very closely and much 
pregraming effort was made to allow for all the boundary 
conditions that Hadge had programmed for. 

The computer program presented here 15 simplified for 
implementation on a VAX computer with one of the standard 


International Mathematical and Statistical Library’s linear 








programing routines. Linear programming 1S used because it 
is fast, inexpensive and easily adaptable to many 
geometries. I Rave allow for only those boundary conditions 
that are frequently encountered in ship structures, that is 
clamped or simple supported. The capability to handle point 
or linear line loads on any grid member has been introduced. 
The ability to handle a uniform load over the entire 
grillage is also important and should be developed. The 
most important aspect ,however, 15 the capability of adding 
additional constraints to the linear programing routines 
that will allow for failures at other than a node. An 
iterative approach to obtain a final solution to this 
general loading scheme is employed. One of the many types 
of problems that this approach is well suited for is 
redundancy study. In developing the capability to introduce 
local failures and the ability to handle large local loads 
the program finds particular suitable use in analyzing ice 
loaded shios structures or slow collisions. I will 
demonstrate its applicability for this by discussing the 
redundancy of a grillage of a ship’s side that 1s loaded as 


one might expect to encounter in ice loaded ship structures. 
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CHAPTER 1 


MODELING CONSIDERATIONS 





The problem of transversely loaded grids covered with 
plating 1S a complex one that, except for very simple cases 
does not have a closed solution. There are several aspects 
of this problem that can be easily simplified without 
Significantly effecting the end results. In ship structures 
it is the goal to obtain a reasonable bound for the ultimate 
strength of a given structure in order to minimize its 
weight and in doing so maximize its efficiency. Tt 1S well 
known that the thin skin plating adds little to the 
transverse strength of a grillage and can easily be included 
as an effective width calculation in determining the fully 
Plastic bending moment of the individual members of the 
Grillage. As such, this paper will deal entirely with the 
strength of the grid. When a grid 1s deflected. twisting 
moments are introduced in grid members. For small 
deformations, these twisting deformations are relatively 
small and will be iqnored. It is also assumed that shear 
force does not contribute to the breakdawn at any plastic 
hinge. For analysis, all grid members are assumed to be 
perfectly elastic and undergoing only small deformations. 
The program will identify a lower bound for ultimate 
collapse and give a possible failure mode. 

The problem reduces to a limit analysis of grillages 


loaded with a variety of transverse loads. In order to 





apply linear programming techniques each load must be 
reduced to corresponding point loads at nodes. This 
reduction does not adversely affect the outcome since it is 
at these intersections that the interactions occur. The 
behavior of the member itself is determined later using the 
moments calculated at each node as the end moments for the 
Given member. Once the load reduction is complete the lower 
bound plastic load limit is found using a standard linear 
pragramming routine. 

As a simple example consider a simple supported nine 
node grid consisting of six members arranged sysmmetrically 
as shown in figure (1). Applying a concentrated load F at 
nede number 3, the center node, will produce a failure mode 
as shawn. The circles indicate fully developed plastic 
hinges. Twisting daes not Jecite in this example because of 
symmetry. This is not the general case. The effect of 
terque in normally encountered structures is relatively 
small when campared to bending. and in ignoring twisting 
moments results will be conservative. Tait (1) has shown 
that for typical I sections the error is less than 0.1 
percent for the types of grids generally considered. For 
long narrow less symmetrical grids the error may be as large 
as 3 percent. 

Determining the work done in creating the deformation 
pattern shown in figure 2., which is the only symmetrical 


failure mode possible far this geometry, the lower plastic 
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Figure 2. A nine node qrid collapse mode 
and deformation profiles 











limit for collapse can be determined. If the angle of 
rotation is ® then the transverse load is Gisplaced by @L 
where Lois the length of each member. And, the external 
work 1S FOL. Each hinge shown undergoes a rotation of @ so 
the total work performed is 4Mo@. By setting the two values 
equal to each other the total collapse load can be 


calculated. 


The program GRIDS presented here produces identical 
solutions. The results of other more complex examples 
completed by Hodge (2) and Tait (1) Rave also been 
duplicated and the results have been identical. GRIDS 1s a 
FORTRAN computer program that uses plastic methods of 
analysis tae determine the response of a grillage subjected 
to a general scheme of transverse loading. The program can 
handle linear line loads on a member, point loads at any 
location an the grillage or a distributed load over the 
entire grillage. This program can be used to find the 
Olastic load limit of any grillage. 

GRIDS is a relatively simple program that employs a 
linear programing technique to solve a system of equations 
that are constrained by the limit of a fully developed 


plastic mament in a member. A program developed by Tait and 





Hodge (1) produced a similar solution but was limited in the 
type of loading that it was able to handle. It allowed only 
point loads at beam intersections. Abbott (2) modified 
Tait’s and Hodge’s program and introduced a Simple geometric 
method of reducing patch loads to nodal loads (loads at 
intersections) and then applied Tait’s and Hodge’s method 
directly to arrive at a solution. GRIDS uses the same 
technique combined with the specific loading scheme to check 
individual grillage members for potential local failures and 
then employ an iterative process in narrowing the bound and 
failure mode based on the points in the grillage not on the 


nodes that are most iikely to fail. 


Pregram Input and Prablem Appraach 


The input data used by GRIDS consists of the specific 
gecmetry. number of beam intersections, beam yield moments, 
lengths of the beams, boundary conditions, and specific data 
for all the loads. The program is limited to one point load 
per node, ane additianal point load per member and ane 
linear line load per member. 

Once all the loading data is entered, equations of 
equilibrium are developed for each member. In reducing 
these equilibrium equations for each member the assumptian 
is made that no local failure of the member has oaccurs. 

This is verified with the output from the first iteration of 
GRIDS. If a particular member is faund that has exceeded 


its yield moment then the moment at the location of this 
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occurrence is constrained. An additional equation 1s added 
to the system for each such member. A new solution 1s 
determined which includes these new constraints. Figure (3) 
shows an equilibrium diagram for a node in a general loading 
condition. 

In order to calculate a solution the following approach 
is followed. The shear forces at the ends of each member 15 
solved for and a moment balance performed. For example. 


taking member (1) we have: 
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The total shear at a node will be the sum of the shears of 
the members that meet at that node. The only unknowns in 
these equations are the end moments. Combining the load 
terms and end shears into a single load for each node vat a2 
solve for the vertical and horizontal moments at each node. 
Referring ta figure (2) we have from equilibrium of 


force at node 1j 3 


ae = V2 + VE + V6 + V7 (3) 


and from equilibrium of moments across each member we have 


VZ = V4 = Cineeieed (4) 


a 


Similar expressions are obtained for all other members. 


Substituting equation (4) into equation (5) we have 


ee 41 at est ede ct Seale nF 
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Mimeitiy = lj 1j 
oe 
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which satisfies equilibrium of node 1j- Similar equations 


are developed for each node. 


The pragram GRIDS constrains each moment to satisfy 





—-Ma < Me. < Ma (6) 


me = 4. < Mo’ (7) 


where Mo 1s the yield moment of a particular member. The 
yield moments of symmetrical beams can easily be calculated 
from the geometry and material properties of the beam. 

Now GRIDS applies a linear programming routine toa the 
system to find a solution. To do this the pragram will 
maximize the load factor PP, where PP is the multiplier of 
each of the applied loads required to induce collapse. In 
MAXIMIZING PP aA Maximum number of maments La and oi; will 
reach the imposed constraint + Mo. 

The linear programming problem seeks to minimize an 
objective function that is a linear function af unknowns, 
subject toa constraints. These canstraints consist of linear 


equalities and inequalities. In the standard form the 


inequalities can be expressed as equalities and the problem 


hecomes: 
minimize £¢(X1,.X%2.2.2222,5A4N) = CX, as C645 ee Mera e CX. (8) 
Neeret@ a..X, + a,~Xr%, © weve te = (9) 
subject toa Aiati + Arak, se Ain*n Dy 9? 
An *y + A554, a vere ots an a = Db. 
anit aE aos + isfactarea ee ae ae = b 


The coefficients of the constraint equations and of the 








objective function are known values and the X'S are 
unknowns. 

Since our problem is to maximize PP 1t should be noted 
that minimizing f (X1,X%2,X%2..22625XN) 15 equivalent to 
maximizing —-f(X1.X2,.K% Sane 2 ee XN). A detailed description of 
linear programing can be found in reference (4). 


The objective function of our problem becomes: 


en. + CX 


11 Fis 


ce ee Oe? ee aed (10) 
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where the X’s are given by 
Xx = M. . +M*M cit) 
n ij Oo 

so that the inequalities (&) and (7) become 
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CHAPTER 2 
PROGRAM OPERATION 

GRIDS is a Fortran program written in a modular fashion 
that makes maximum use of separate subroutines and common 
memory. This approach makes the program much easier to 
understand and gives it the flexibility to be easily 
changed. 

The program has some specific limitations. All 
horizontal members and their maximum yield moments must be 
the same. All vertical members and their maximum yield 
moments must be the same. However, vertical and horizontal 
members need not be the same. Boundary conditions may only 
be clamped or simple supported. But, the boundary condition 
of each side 1s independent of the others. Point loads may 
be applied at any point of on the grillage, but there must 
be cnly one point lead per member with one additional point 
load allowed at each node. One linear line load per member 
1s allowed. All loads must be applied perpendicular and 
transverse ta the grillage. 

The main program GRIDS assimilates major variables and 
input and sequences the subroutine operation. There is only 
one set of significant calculations performed in the main 
program. The calculation consists of determining the 
dimensions of the arrays that the IMSL linear programming 
routine ZXZLP uses and could not be done in the subroutine 


that actually calls ZXSLP. 





Each of the program’s functions is controlled by an 
independent subroutine. The controlling program GRIDS takes 
the user input for the grillage to be studied and assigns 
variables for the boundary conditions and geometry of the 
structure. subroutine INTPL initializes the load array. 
MINT uses the specific boundary conditions and initializes 
the output load arrays HMNT(i1) and VMNT(i). The index i is 
for the 1th node. These arrays are filled with ones and 
zeros. All nodes that are known to have a zero moment 
because of the boundary conditions are assigned a zero. all 
other unknown moments are assigned a one. These arrays are 
used to determine the total number of unknowns and are used 
as multipliers in Subroutine CONST to ensure that the 
boundary conditions are handled properly (ie. zero moments 
at simple supports). Later they are used to store the nodal 
moments calculated in the linear programming routine. 

Subroutine ENTRL handles assiqning all loads to the 
preper arrays and decomposes them into nodal loads for use. 
Fach load must be stored so that individual beam moment 
distributions can later be calculated by Subroutine BEAMCK. 

Once the loads are entered the user decides if yield 
moments are allowed to form in the members away from the 
nodes. If he chooses to allow additional plastic hinges 
then the subroutine CONST is called and the additional 
constraint equations are written. These equations may be 


generated automatically in which case they are based on the 





grillage scheme and an initial check of each member. They 
can also be generated one at a time manually. In this case 
the user inputs which member is to have the additional 
constraint and the location on the member that is expected 
to develop the plastic hinge. In addition, if the manual 
mode is chosen, the constraint equation can be written 
normally as a likely failure point or it can be written as 
if the point has already failed due to some very large load 
(in this case it becomes an equality constraint). These 
possibilities are the heart of the program and are discussed 
in detail in the body of this paper. 

Subroutine EONS is used to calculate the nodal 
equilibrium equations that are sent to IMSL’s ZX=LP. 
Ssubrautine Fixer transforms all of the equality and 
inequality canstraints into arrays that ZXSLP can use and 
calis ZX2LP. Tt then translates the results into horizontal 
and vertical maments at each node. 

Once this is complete the moments at each node and the 
Specific loading scheme are used by RBEAMCK to determine the 
maximum moment and its location on every loaded member. If 
automatic constraints was chosen‘at the beginning of 
execution then the program will automatically constrain any 
moment found to be greater than the yield mament and iterate 
antil a solution is found where no moment exceeds the yield 
moment or until the maximum number of iterations 1s reached. 


Subroutine OUTPUT writes all grid information to a data 
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file named MOMENTS.DAT. All inputs are summarized and a 
ZXSLP error number is provided. Output includes all members 
ind moments, nodal loads (from decomposed general loads) and 
the maximum moment and location for each loaded member. 


Cc 


Figure 3 shows the macroscopic program flow. 


SUBROUTINES 
MINT AND INTPL 
MINT initializes the nodal bending moment arrays. Any 


moment that is known to be zero because of boundary 
conditions are set to zero and all others. which are 
unknowns, are set to one. The master program GRIDS then 
Calculates the number of unknowns in the problem by counting 
the non-zero elements in the array. These ones and zeros 
are used as multipliers to cause the coefficients in the 
equilibrium equations calculated in CONST to go to zero toa 
comply with boundary conditions. 

INTPL sets the starting value of each nodal load toa 


.2FroQ. 


ENTERL 

ENTERL is the subroutine that makes load variable 
assignments. Its cperation is straight forward and requires 
little explanation. Some error traps are built in to 
ensure all loads are in fact valid from the standpoint of 
location. There is no trap however to ensure that the 


restrictions on the numbers of various loads are observed. 











If these restrictions are not observed the calculated 


Arrays 


PTLOAD (1) 
II di) 
JJ (1) 
{iL DIND (1) 


DIVAL (1) 


DILOC (1) 


LD2ND (13) 


D2AVL (1) 


D2RBVL (i> 


i DSND (1) 


DEVAL (1) 
VMNT (1 > 
HMNT (1) 


RHEMNT (3 > 


RVEMNT (3) 


PROGRAM VARIABLES 


The total load for node (1) 

The Y coordinate of node (1) 

The X coordinate of node (i) 

The reference node number of point load (1) 
The magnitude of point (41), a negative 
value indicates the point load is on 
vertical member. a positive value 

indicates the load 1s on the horizontal 


member 


The distance of point load (1) from the 
reference node 


The reference node number of line load (1) 
The magnitude of line load (1) at the 
reference node. a positive value 

indicates the load 1s on the horizontal 
member 

The magnitude of line load (41) at the node 
opposite to the reference node, a positive 
Value indicates the load 1s on the Horizontal 


member 


The reference node number of distributed 
load (1) 


The magnitude ef the distributed load 
The value of the vertical moment at node (1) 
The value of the horizontal moment at node (1) 


Maximum moment on horizontal member of 
reference node (41) 


Maximum moment on vertical member of 
reference node (41) 


Table 1 Program Variables 





RVELOC (1) 


RBHELOC (2) 


ACNST (1) 


NCNST (1) 


FOG(i.j)\ 
Rei. j) \ 
Bi (i) / 
B2 (i? y 


Variables 


NODESH 


NODESV 


NODTOT 


VLEN 


HLEN 


ERY 


LNLDN 


tDIS 


ITOP 


ILFT 


Table i 


Location of max moment on vertical member of 
reference node (41) 


Location of max moment on horizontal 
member of reference node (1) 


Location of manual constraint for 
reference node (1) 


Reference node of manual constraint (1) 


These are dummy arrays for manipulating data 
for use by the linear programming subroutine 


Number of beam intersections in the horizontal 


direction (nodes) 


Number of beam intersections in the vertical 
direction (nodes) 


Total number of beam intersections. all edges 


have beams along them 


Length of the vertical members. distance 


between nodes 


Length of the horizontal members. distance 


between nodes 
Number of point loads added to the aqrillage 


Number of linear line loads added to the 
grillage 


Number of distributed loads added to the 
grillage 


Top edge boundary condition 
Left edge boundary condition 


Program Arrays and Variables 
(continued) 





IRGT 


IRITM 


PP 


IER 


TEQNM 


NODEE® 


HORZMO 


VERTMQ 


Right edge boundary condition 

Rottom edge boundary condition 

Grid load factor 

Linear programming routine error number 
Number of non-zero moments in the grid 


Number of nodal equations required to 
obtain a solution 


The fully plastic moment of the horizontal 
Grid members 


The fully plastic moment of the vertical 
Grid members 


Variables used throughout the subroutines are 


dummy variables used in calculations or for keeping track of 
vector sizes and locations. 


Table 1 Program Arrays and Variables 


(continued) 





maximum moment and location on each member will not be 
correct. 

There are several arrays for each load type. The total 
number Of each type of load is kept track of with load 
counters. LPT is the total number of point loads entered 
and LNLDN is the number of line loads. LDIND and LD2ND 
store the reference node numbers of each point load and line 
load respectively. D1VAL and D1ILOC store the magnitude and 
location (distance from the reference node) of each point 
load. D2AVL and D2BVL contain the end point magnitudes of 
each line load. 

Table (2) Gives a summary of load array functions. 

All array elements are indexed by load number for each 
specific type of load. Fach load type has an array that 
contains the reference nodes associated with the load 
numbers. The sign of the load magnitude determines which 
member, either horizontal for positive or vertical for 
negative, the load is associated with. The reference node 
is always the node to the left or above the loaded member. 


Figure 6& shows the program flow for entering loads. 


EQNS 





Subroutine EQNS writes the nodal equations of the 
Qrillage and writes them to an array ina format that 
subroutine FIXER uses later to assimilate all of the 


required data into the proper format. EQNS caiculates the 





Point Loads LDIND (1) 
DIVAL (i)” 
DILGE Cs) 

tine Loads L_DOND i) 


DZAVL (i)” 


_.# 
D2RVL (1) 


odal Loads PTLOAD (1) 


Contains the reference 
node number for point load 
number 1 


Contains the magnitude of 
load number i 


Contains the location of 
(distance from the 
reference node) load 
number 1 


Contains the reference 
node number for point load 
number i 


Contains the magnitude of 
load number 1 at the 
reference node 


Contains the magnitude of 
load number 1 at the 
opposite node 


Store the total nodal 

load (decomposed general 
loads) that the equilibrium 
equations are written from 


These values are set to the negative to indicate the ioad 
is on the vertical member associated with the reference node 
indicate the load is on th 


and are left pasitive to 


horizontal member. 


Table 2. 


Load Array Function 









CHOOSE TYPE OF LOAD 
TO ADD 

POINT OF LINEAR LINE 

LOADS 





RETURN 















INPUT: REF. NODE (NDNUM) 
HORZ OR VERT MEMBER (DIR) 
MAGNITUDE (VAL OR AVAL OR RVAL) 
OPPOSITE NODE MAGNITUDE OR LOCATION 
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—— 
we | 


IF NOT 
OK 


DECOMPQSE TO NODAL LQADS AND 


ADD TO FTLOAD ARRAY 
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Figure 6 Program Flow for Entering Loads (ENTERL) 





coefficients of the equilibrium equations for each member. 
Since the unknowns are the end moments, a single linear 
equation that contains each unknown throughout the entire 
Grillage is generated for each member. The result is a 
vector with twice the number of nodes plus one elements. 
since each member has only two ends, all of the coefficients 
except two turn out to be zero. Only the two coefficients 
associated with the nedes that bound the member will be non- 
zero. The additional unknown is the load factor which is 
not censtrained but will be maximized. The left half of 
this array represents all of the horizontal moments and 
members while the right half represents the vertical moments 
and members. Figure (7) shows the construction of the 
yvecter matrix for all equations written by both subroutines 
EQNS and CONST. 

In Figure (7) the first example equation represents the 
equilibrium equation to the member located between the node 
in the third and forth columns ar node numbers three and 
four, The last column represents the load factor and the B 
vector represents the right hand side of the equation. 

Since the non-zero elements are located in the left half of 
the matrix we knew that this equation represents a 
horizontal member. The second equation has non-zero 


elements in the right half of the matrix so 1t must be for a 





a eee el AZ Matrix --<-------------- > 
Jota €--total # nodes-->:<--total # nodes-->: 
+ 1 2 Se. Oh. eee isl: oc Bels Geemarc eats et) = 
tt Horizontal Section : Vertical Section : 
n z F 
i Example equations : : 
n % = 
oO OF Oe ee Omen Ors Oro 1O OLONe GCG. O20 "5 
Ww Oo OF OO Are OOO 3520 O° 0 2 OLe7O 2 © & 
n = = 
Ss = rs 
: ; = 
\il = 3 
Figure /. 


vertical member. Since all of the nodes of the grillage are 
numbered from left to right we know that all of the 
horizontal members will be represented by coefficients in 
adjacent elements. Likewise. vertical members will be 
represented by elements in the right side of the matrix that 
are separated by the number of horizontal node in the 


Qrillage. Fiqure 8 shows the flow diagram for EQNS. 


CONST and CONST2 

Subroutine CONST constructs additional constraint 
equations in order to refine and clarify the problem 
solution. Specified members are selected for these 
additional constraints. Coefficients of the equilibrium 
nodal equations for each of these specified members in the 


grillage that has unknown end moments are then calculated. 
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; RETURN 


Figure 8 EQNS Subroutine Flow Diagram 








It fas Automatic or manual options. If automatic 1S 
selected the program calculates the maximum moment found in 
a member as the sum of the linear combination of the 
member*’s end moments as calculated from the initial nodal 
loading and the particular loading of the member. The 
coefficients of the nodal equations involving that member 
are then calculated using the location of the maximum moment 
found during this initial look at the member. In manual, the 
location of the moment is chosen by the user. In manual the 
user also chooses whether or not the extra constraint 1s to 
Have a moment already developed (1e prior damage) or it 15 
to be handled as a limit. 

CONST2 operates exactly as subroutine EQNS ta write an 
Array of Coefficients to be used by subroutine FIXER. 


Figure 9 shows the flow diagram for CONST and CONST2. 


FIXER 

Subroutine FIXER gathers the known data and converts it 
into a form that can be used by the IMSL linear programming 
routine, ZXELP. It writes an identity matrix that acts as 
the list of inequality constraints and it uses the arrays 
produced in EQNS and CONST2 to write the equality 
constraints and additional constraints. All of these are 
combined in a single matrix that is used by ZXSLP. see 
Appendix B for details of this matrix. ZXSZLP is then 


called. FIXER then converts the solution returned from 


7 





WRITE INEQUALITY CONSTRAINTS 
Xn < Mo 
IN FORM OF IDENTITY MATRIX 


CONVERT NODAL EQUILIBRIUM 
EQUATIONS FROM MATRIA FOG AND BI 
TO CORRECT FORMAT 


CONVERT ADDITIQNAL CONSTRAINT 
S<QUATIONS FROM MATRIX AZ AND B2 
TQ PROPER FORMAT 





CALL LINEAR PROGRAMMING 
ROUTINE ZXSLP 


RETURN 


TRANSLATE RESULTS TO MOMENT FORMAT 


Figure 10 Subroutine FIXER Flow Diagram 








a) 


BEANCK 

Subroutine BEAMCK calculates the maximum moment and 
location for every loaded member. It accomplishes this by 
looking at each node and gathering the loads that are 
referenced to it. If no loads are referenced there will be 
no calculations for that node. It then separates the loads 
into vertical and horizontal members and assigns them to 
dummy variables. Equations for the moment as a function of 
position are available for each possible loading condition 
and position on a member. There are four possibilities. 

The linear line load may be either increasing or decreasing 
GQInNg away from the reference node and the moment must be 
Calculated for either side of a point load. A constant line 
load will fir either type of linear line load equation. 
These equations are derived in Appendix C. 

The program applies the proper equation and calculates 
the moment and location. The results are stored in four 
arrays. BHEMNT and BHKLOC store the moment and location for 
the horizontal member referenced to a node and RBVEMNT and 
RVKLOC store this information for the vertical members. It 
Should be noted that subroutine CONST uses identical logic 
to determine which additional constraint equations to write 
when automatic constraints is selected. Figure 11 shows the 


subroutine flow for BEAMCE. 


OUTPUT 


Subroutine QUTPUT consists almost entirely of write and 





GATHER ALL THE LOADS 
ASSOCIATED WITH A NODE 
AND ASSIGN TQ DUMMY VARIABLES 
IF NO LOADS SKIP THE NODE 


CHOOSE THE CORRECT EQUATION 
FOR THE LOADING CONDITION ENCOUNTERED 


USING THE END MOMENTS CALCULATED 
USING ZXSLP, CALCULATE THE MAXIMUM 
MOMENT AND LOCATION 


STORE THE MAXIMUM MOMENT AND LOCATION 
IN BHEMNT AND BHKLOC (HORZ MEMBER) 


OR BVEMNT AND BYKLOC (VERT MEMBER) 


RETURN 





Figure 11 BEAMCK Subroutine Flow Diagram 





format statements that provide a hard copy of the grid 
Characteristics, loading and the solution. A ZXSLP error 
number iS also provided which can be referenced in Appendix 
R if problems are encountered. A sample output 1s shown in 


Rigure (Lee 


RUNNING THE PROGRAM 

Once GRIDS and its Subroutines. including the IMSL 
routine ZXZLP and ZX20P (ZX20P is used by ZX=3LP), are 
installed on the operating system, the program can be run. 
It may be run interactively using keyboard input or it may 
be run using an input data file. Appendix B shows how to 
create an input data file. 

The following listing is the actual program dialogue as 
seen on the terminal when running Grids. User responses Rave 
been put in square brackets for clarification. The listing 
below is typical but not all inclusive. If automatic 
constraints is selected, the inputs that concern manual 
constraints is not asked for. AIL1l other inputs are similar. 
C$ RUN GRIDSI 
PROGRAM RESTRICTED TO A MAXIMUM GRID OF & X 53 
IE 30 NODES TFOTAL OR LESS 
THIS RESTRICTION IS IMPOSED BECAUSE OF THE MAXIMUM 
SIZE OF THE WORK VECTOR USED BY ZxX3LP IS CONSTRAINED 
RY THE PROGRAMMER. 


SEE PROGRAM MANUAL FOR MORE INFO 


INPUT NUMBER CGF NODES ALQNG THE TGP AND SIDE OF THE GRID 


C4.4] 
INPUT BOUNDARY CONDITIONS FOR THE TOP 





LEFT.RIGHT AND BOTTOM EDGES OF THE GRID 
IN THAT ORDER 

INPUT A 1 FOR SIMPLY SUPPORTED EDGE 
INPUT A 2 FOR A CLAMPED EDGE 

ALL VALUES ARE INTEGERS 


Mtwde 1,13 
INPUT HORIZONTAL MEMBER LENGTH (REAL NUM.) 


Lee 
INPUT VERTICAL MEMBER LENGTH (REAL NUM.) 


Cis. J 
INPUT HORIZONTAL BEAM YIELD MOMENT (REAL) 


C20000, ] 
INPUT VERTICAL BEAM YIELD MOMENT (REAL) 


C50000. 1] 
LQADS ARE ENTERED REFERENCED TO THE NODE 
ABOVE OR TO THE LEFT OF THE LOAD 
ALLCWED LQADING CONFIGURATIONS: 
1. POINT LOADS 
2. LINEAR LINE LOADS ALOQNG A MEMBER 
&. LOAD ON A PLATE ELEMENT 


ENTER THE TYPE OF LOAD TO ADD 

POINT LOQADS 

LINEAR LINE LOAD ALONG A MEMBER 
NO MORE LOADS TQ ADD 


D> PJ ee 


Ei 

THIS ROUTINE ADDS POINT LOADS TO THE GRID. THERE MAY RE 
A MAXIMUM QF ONE POINT LOAD PER NODE PLUS QNE PER MEMBER 
THE REFERENCE NODE IS ABOVE OR TO THE LEFT OF THE LOAD 
ENTER FOUR VALUES TO DEFINE EACH LOAD 


IST VALUE = REFERENCE NODE (INTEGER) 

2ND VALUE = VERTICAL OR HORIZONTAL MEMBER 
1 = HORIZONTAL MEMBER (INTEGER) 
2 = VERTICAL MEMBER CINTEGER) 
Q = ON THE NODE (INTEGER) 

“RD VALUE = LOAD MAGNITUDE (REAL) 

4TH VALUE = DISTANCE FROM THE REFERENCE 


NODE (REAL) (IF THE 2ND VALUE = 9 
THEN THE 4TH MUST = Q) 
ENTER ALL ZEROS WHEN COMPLETED WITH POINT LOADS 


(45.1.7-5,7-5] 
ENTER FOUR VALUES TQ DEFINE EACH LQaD 

1ST VALUE = REFERENCE NODE (INTEGER) 

2ND VALUE = VERTICAL OR HORIZONTAL MEMBER 








1 = HORIZONTAL MEMBER (CINTEGER) 
2 = VERTICAL MEMBER (INTEGER) 
O = ON THE NODE (INTEGER) 

SRD VALUE = LOAD MAGNITUDE (REAL) 

4TH VALUE = DISTANCE FROM THE REFERENCE 


NODE (REAL) (IF THE 2ND VALUE = O 
THEN THE 4TH MUST = 9) 
ENTER ALL ZEROS WHEN COMPLETED WITH POINT LOADS 


[0,0,9.,9.] 
ENTER THE TYPE OF LOAD TO ADD 

POINT LOADS 

LINEAR LINE LOAD ALONG A MEMBER 

NO MORE LOADS TO ADD 


2 PJ 
wou on 


CO] 

THIS SUBROUTINE ALLOWS PLASTIC MOMENTS TO BE 
FORMED AT LOCATIONS OFF OF THE NODES OF THE 
GRILLAGE. THIS IS DONE BY WRITING ADDITIONAL 
CONSTRAINT EQUATIONS FOR THE LINEAR PROGRAMMING 
ROUTINE 

AUTO CONSTRAINTS WILL LOOK AT EACH MEMBER AND 
CONSTRAIN THE MAXIMUM MOMENT TO <= YIELD MOMENT 


MANUAL ALLOWS THE USER TO CHOOSE WHICH MEMBERS 
TO SET CONSTRAINTS FOR AND THEIR LOCATION 


INPUT A 1 FOR AUTOMATIC CONSTRAINTS 
A 2 FOR MANUAL CONSTRAINTS 
A O FOR NO CONSTRAINTS 


C24 

CONSTRAINED MEMBERS MUST HAVE AT LEAST ONE LOAD 
ASSOCIATED WITH THEM 

INPUT THE REFERENCE NODE FOR THE ADDITIQNAL 
CONSTRAINT. AND A 1 FOR HORIZONTAL MEMBER 

OR A 2 FOR VERTICAL MEMBER (INTEGERS) 

ENTER ZEROS FOR NCQ CONSTRAINTS 


45,1] 
INPUT THE DISTANCE FROM THE REFERENCE NODE 
TO LOCATE THE CONSTRAINT (REAL) 


(7.9) 

ENTER A 1 TO INDICATE THAT THE MEMBER HAS A 
FULLY DEVELOPED YIELD MOMENT ALREADY DEVELOPED 
NOTE: THIS TYPE MUST BE ENTERED LAST 

ENTER A O IF THE HINGE IS NOT ALREADY FORMED 


{ INTEGERS) 


COJ 





At this point the program will run. All of the output 15 
written to a data file that can either be written to the 


terminal or printed on the system hardcopy device. 





REF 


RETWEEN NODES 


(2) oS [6p [8 ), 


-HHORZ BEAM PLASTIC BENDING MOMENT 


GRID CHARACTERISTICS AND LOADING 


HORIZONTAL LENGTH 15.90 
VERTICAL LENGTH 15.900 


NUMBER OF NODES HORIZONTAL 4 


NUMBER OF NODES VERTICAL 


4 


BOUNDARY CONDITIONS 


SIMPLE SUPPORTED 
SIMPLE SUPPORTED 
SIMPLE SUPPORTED 
SIMPLE SUPPORTED 


TOP 
ROTTOM 
BEET 
RIGHT 


O. 10E+04 


O. 10E+04 


VERT BEAM PLASTIC RENDING MOMENT 


ZXZLP ERROR NUMBER = QO 


POINT LOADS ADDED TO THE GRID 


NODE DIR FM NODE DISTANCE FROM NODE MAGNITUDE LOA 


NONE USED — 


LINE LOADS 


REF NODE MAG OTHER NODE MAG LINE LOAD NUI 


AND v4 1.00 ze OO 1 
AND 10 1.00 22 OO 2 
AND 11 22Q0 1.Q0 = 
AND ii 2200 1.00 + 
NODE LOADS USED IN CALCULATIONS 


NODE NUMBER NODE LOAD 
0.00 
QO. QO 
QO. QQ 


Led Boe 


Figure 12 Sample Output 





VONCGUSD 


10 
11 
ie 
13 
14 
15 
16 


bh) KJ 
290 
o 


2 DOs) OsD Die oe | 
oo oO OS 2 2 


NJ BJ 
Se OOS F9OUN DONS 


Oo Ca & 


ADDITIONAL CONSTRAINTS USED 


RETWEEN NODES 


QO AND 


*e*#*%* ITERATION NUMBER IS 


TYPE USED 
Q NONE 


—stss 


NODAL MOMENTS 


NODE NUMBER 


THE GRID LOAD FACTOR = 


SOOnNOM AW he 


HORIZ MNT 
QO. Q0E+00 
QO. OOE+O00 
0. Q0E+00 
O. QOOE+00 
O. OOE+00 
O.10E+04 
Q. 10E+04 
0. O0E+00 
9. O0E +00 
QO. 10E+04 
QO. 1OE+O04 
O, OOE+00 
Q.Q0E+00 
O. QQ0E+00 
2. QOE+00 
O.Q0E+00 


> 


ss. 


Je 7k 


VERT MNT 


Q.00E+00 
O. QOE+00 
Q. 00E+00 
OQ. Q0E+00 
0, Q0E+00 
OQ. 86E+0% 
Q.10E+04 
Q.O0E+00 
Q. O0E+00 
QO. 10E+04 
OQ. 86E+02 
Q, OOE+00 
Q. OOE+00 
- NOE+00 
OE+00 
OE+00 


OO 
o> 


MAXIMUM MOMENTS AND LOCATION ON LOADED MEMBERS 


BETWEEN NODES MAGNITUDE 
HORIZONTAL MEMBERS 

& AND 7 1242, 

10 AND 11 1242. 
VERTICAL MEMBERS 

6 AND 10 1179. 

7 AND 11 1179. 


**% INDICATES 


Figure 12 Sample Output 


> YIELD MOMENT 


LOCATION 


(continued) 





CHAPTER 3 
DISCUSSION and EXAMPLE PROBLEM 

Work done by ABBOTT (2) resulted in a program capable 
of estimating grid load factors based only on nodal loads. 
These estimates turn out to be overly conservative. For any 
symmetrically loaded grillage one can expect the failure mode 
to also be symmetric. In applying linear programming to find 
a failure mode (or any technique that assumes nodal loading)? 
a load factor and a pattern of plastic yield hinges at the 
Grid intersections is developed. It should be noted that 
the solution found will not necessarily be unique. 

In maximizing the number of plastic hinges for any 
given system there is bound to be adjacent nodes with fully 
developed plastic hinges. This will be true in all cases 
except the nine node grid. If a member has two fully 
developed yield moments of the same sign at each end, then 
the moment developed within the span must be greater than at 
the ends for any cansistent loading on the member. The 
implication is clear. For any consistent loading scheme 
there must be additional yield moments formed at locations 
ether that at the nodes. In comparing the results of the 
program GRIDS to the results presented by Hodge (1) or 
Abbott (2) it can be seen that a program that allows for 
yield moments to be formed with in a span returns a much 
more acctrate result. 


In taking a simple symmetric case and analyzing it for 





various conditians the use of GRIDS can be demonstrated. 

The program’s limits and a straight forward verification also 
result from this exercise. A nine node grillage was 

analyzed using a point load to compare relative load 

carrying capability and collapse modes for various boundary 
conditions and load conditions. 

The first geometry considered is the nine node grillage 
with clamped boundaries. The load was first applied at the 
node and then gradually moved toward the edge along one 
member. First. if this grillage 1s analyzed in the same way 
as the simple supported grid in chapter two was, the grid 
load factor is faund to be 8 A yield moment forms at the 
intersectian of the members in both the horizontal and 
vertical directions. With the edges clamped hinges will 
also farm at the boundary end of each beam. Setting the 
work dane in deflecting the node to the sum of the work done 
in forming each of the plastic hinges will give the upper 
Plastic limit or collapse load of the arid. The following 


equalities are obtained. Here F is the collapse load. 


8 Moa (14) 


aT 
- 
D 

It 


F = 8Mo/~L (15) 








Figure 12 Collapse Mode With Plastic Hinge at Midspan 





Now, 1f the point load 15 moved to midspan. we can 
calculate a new collapse load and allow for a plastic hinge 
to form at the point on the member that will see the maximum 
moment. Figure 14 shows the collapse mode and load 
configuration for a point load located half way between the 
node and the edge of the left center horizontal member. 
Because of the symmetry involved in this case 1t 15 
relatively simple to calculate the collapse load by hand. 
Summing the total work done in deforming the grid we Rave 
the force applied moving through a distance of QL/e and the 


center node with its effective force displaced a distance of 


Q.. 
1 


F/2 @.L + FO, L/2 (16) 


{. 
| 
boob 


And summing the work required to form the plastic moments 


Gives: 


zMoe,. + 2M09 + 6&Moe (17) 


Setting (3) and (4) equal to each other and realizing that 


AX is the sum of 8, and 8, gives: 


i = 


FL(26, + @_) = 8Mo (20, + 6.) (18) 


F = 8Mo0/L Cir} 





The fact that this result is the same as found in (15) 
above is only coincidental and is due to the symmetry 
involved. Similar calculations can be performed for any 
position along a member. A tabular comparison of the 
results of the program Grids and the results of a program 
that allows only nodal loading as a function of location of 
the point load is presented in Table 3. For the application 
of a point load and no intermediate yield moment, it is 
readily apparent that the farther the load moves from the 
node the more error is introduced in computing the collapse 
lead. This occurs because the yield moment that must form 
at the point load will be formed with much less load as the 


load moves away from the node. 


CLAMPED NINE NODE LOAD FACTORS 


LOCATION « LOAD FACTOR WITH LOAD FACTOR WITHOUT 
IN 4 OF PLASTIC MOMENTS PLASTIC MOMENTS 
LENGTH IN THE MEMBER IN THE MEMBER 
10 76/.68 888.89 
sO 747.29 1142.86 
39 Boo. 1600. 

70 9352.58 2666.67 
70 eee Boo0Q. 


* measured from the center node 


Table = Load Factor for Various Load Locations 


For Plastic Moments in a Member and For 
no Flastic Moments Except at nodes 





If the load is normalized to a total nodal load of one 
for each position along the member, the relative load 
carrying capability of the grillage can be plotted as a 
ratio of the load factor at a point divided by the load 
factor for nodal loads only versus position along a member. 
Table 4 lists the data for the nine node example and Figure 


14 shows the result graphically. 





YIELD HINGE IN MEMRER NO YIELD HINGE IN MEMBER 

Ee LOCATION (7) oe LOCATION (7%) 

alge FROM CENTER il FROM CENTER 
NODE NODE 

1 2) 1 O 

0.8636 10 1 10 

0.6328 20 1 20 

0.35 wO 1 DO 

0,3571 790 1 79 

0.2778 90 i 90 


TABLE 4 Normalized Load Factor VS Position of Hinge 


1 — 
.8 - 
PP 
PP ~o - 
n 
4 paw 


10 2O 2O 4D 20 60 7Q BQ 90 


Location (per cent of length) 


Figure 14 Normalized Collapse Load VS Position of Hinge 





The second type of problem tested was a grillage that 15 
a more realistic model of a ship’s side. A6& X 5 grillage 
with a line load was analyzed. A line load is a reasonable 
approximation of a ridge of ice along a ship’s side. The 
kine lead was placed along a series of horizontal members 
and the load factor and failure mode were calculated. The 
load was moved vertically downward and the calculations were 
repeated for different positions along the vertical beam. 
Note that once the line load is off of the horizontal 
member, it acts like a series of point loads along a line 
across each of the vertical members. Figure 15 shows the 
loading and a typical failure mode encountered. Table 53 
shows the camparison of the load factors given by grids to 
those that are the result of a program that allows no 
intermediate hinges to form. It is interesting to note that 
the load factor remains more nearly constant when 
intermediate yield moments are allowed to form. 

The final problem to be considered is one that was 
considered by Abbott (2). Here we must convert the 
distributive load into line loads on the members. The 
prablem is an example that illustrates how GRIDS might be 
used in considering the strength of a ship’s side that 1s 
anticipated to operate in a particular ice environment. For 
this problem consider that the ship 1s to operate in a 
uniform ice sheet five feet thick. Further assume that the 


ice has a ltacal crushing strength of 300 psi and that 





LOCATION (%) WITH INTERMEDIATE NO INTERMEDIATE 

PLASTIC MOMENTS PLASTIC MOMENTS 
Oo 384.00 400.00 
10 380.36 412.79 
30 nace as 444.44 
50 337.78 444.44 
79 330.78 292.16 
90 330.8% 350.88 


Table S Load Factor for 5 xX 6 vS Position of Line Load 





®=POSITIVE YIELD MOMENT 
@ =NEGATIVE 


Figure 15 Failure Mode for Line Load on Horizontal Member 
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with a Horizontal Line Load Between 
Horizontal Members With Clamped Boundaries 





crushing 1s the predominant failure made. Analyze a ship 


which has the following characteristics. 


Frame Spacing zi inches 
Rulkhead Spacing 7 feet 

Distance Retween Decks 1S feet 

Number of longitudinals z (evenly spaced) 
Stringer Modulus a a 

Frame Modulus -8 in® 

Yield strength YJOOOO psi 


Figure 17 shows the configuration for this problem. 
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Figure 17 Ice Loaded Ship’s Side 





Converting the distributive ice loads on this grillage 
we get lineloads around each plate element. Each plate has a 
total pressure of 3.78 X fom pounds on it. Assuming this 
reduces to a uniform line load around the perimeter of each 
plate results ina line load of 583.32 pounds per inch. 
Entering the data and selecting automatic constraints 
results gives a Load Factor of 9.175. With this load factor 
we know that this grillage can support less than twenty 
percent of the anticipated load. The structure must be 
strengthened. Strengthening can be accomplished in several 
ways. Increasing plate thickness, frame of longitudinal 
member sizes or using stronger materials are all viable 
options that will result in an increased maximum yield 
moment of each member. Increasing the yield moment of each 
member results in a grillage that will support the applied 
load. 

Additional insight can ge gained by varying the strength 
of either the horizontal or vertical members alone and 
ebserve the structural behavior. With the same load applied 
as in the problem above, the horizontal yield moment was set 
at a value of 100,900 and the vertical bending moment was 
varied through a range of values of 20,000 to 90,000. Al}l 
of these values give the same Load Factor of 9.18. This 
suggests that for the applied loading the grillage 1s much 
weaker in the horizantal direction, which makes sense when 


the geometry is considered. In looking at the failure mode 





provided by GRIDS we observe that plastic hinges are 
indicated for all boundary nodes. We also know that there 
are a number of yield moments given that do not actually 
occur, due to the degree of indeterminancy and linear 
programming methods (Chapter 4 contains additional 
discussion on this point). We can conclude from the above 
that the moments in the vertical members are the ones that 
are not actually formed. 

Now, if we hold the yield moment of the vertical 
member constant and vary the yield moment of the horizontal 
member we can observe the grillage behavior as a function of 
the horizontal member size. Figure 16 shows how the load 
factor varies as the horizontal member yield moment for this 


example. 
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Figure 198 Load Factor Versus Horizontal Member strength 
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CHAPTER 4 
SUMMARY OF RESULTS 

As illustrated in Chapter 3 the results that Grids 
gives are significantly more accurate than for a program that 
allows only nodal loading. In choosing to use linear 
programming to solve a problem with intermediate yield 
hinges, an iterative approach is also required. If one 
were to write constraint equations for forming the plastic 
hinges within a member directly, as a function of a general 
combination of point and line loads, the resulting equations 
would be non-linear. In order to use linear programming and 
Still allow plastic hinges in a member, the end moments must 
be used. 50, the problem is first solved using only nodal 
loads to calculate the end moments. Constraint equations 
are then written that are a linear combination of loading 
and end moments. The assumption is that the final end 
moments will not change the basic shape of the moment 
distributian. The end moment only changes the slope of the 
line that the moment distribution lies on. Fram this one 
can see that if the end moments. as given by GRIDS, for a 
particular member are equal then the location of the maximum 
moment is exact. Hence the position of the yield moment in 
that member 1s correct. 

For point loads only the position of the failure is 
known and the program solution 1s accurate. In the case 


of line loads GRIDS gives a "linearized" approximation. 





Calculating the approximate error can be accomplished by 
comparing the location of the maximum moment based on the 
moment distribution and member end moments to the position 
that the additional constraint was written for. Doing this 
assumes that the yield moment that forms at the location of 
the constraint does not change the moment distribution. For 
a constant line load, the location error in percent of the 
length of the member is approximately one half of the 
percent difference between end moments. In other words, if 
the end moments are within twenty percent of each other then 
the plastic yield moment will be within ten percent of the 
locatian given by GRIDS. Further more. it will be shifted 
toward the end of the member with the larger moment. 

The location that GRIDS calculates for the maximum will 
be exact for point loads and constant line loads ona 
sysmmetrical grid with symmetrical loading. In non- 
symmetrical grids with multiple load types the locations of 
the maximum moments are approximate but very close to 
actual. This occtirs because of the way in which linear 
programming calculates the end moments of each member. If 
the end moments of a member are unchanged between the first 
and second iterations as labeled in the output. then the 
location of the maximum moment in that member calculated by 
grids is exact. If one of the end moments of the member is 
reduced then the effect is to lower that end of the moment 


curve until the maximum moment as canstrained 1s found. The 





curve keeps its original shape however. The error 
introduced is small as long as the end point magnitudes of a 
member”’s linear line loads do not differ more than fifty 
percent. 

GRIDS offers both automatic and manual modes for 
developing the additional constraint equations that actually 
add the plastic yield moments in a member. This gives added 
flexibility in evaluating a particular grillage. By 
choosing the manual mode and choosing equality constraints a 
failure can be force to form at a particular location. By 
forcing a failure, or by having an already formed plastic 
hinge in a member, the load carrying capability of a 
partially failed grid can be examined. This technique may 
be desirable to evaluate a grillage@ that has undergone a 
lecal impact on a member so that a hinge may be formed at a 
location first that might otherwise not reach its yield 
moment until much higher general loads are applied. if 3 
however, a plastic hinge is placed ina member with no loads 
associated with it the output is meaningless and a ZXSLP 
linear programming error will occur. 

As discussed previously, the failure mode provided by 
GRIDS is not unique. All of the systems encountered contain 
more unknowns than equations and therefore cannot be solved 
directly. The way in which this linear system is solved and 
Maximized causes selected inequality constraints to be 


evaluated as equality constraints. The number used in this 





fashion is equal to the degree of indeterminancy of the 
system. This means that not all of the locations that are 
listed in the output may actually reach the yield moment. 
But, the output provided is a maximized and feasible 
solution. In other words, there may be fewer yield moments 
farmed in an actual failure but the total work done is 
correct. This could be accounted for by different 
magnitudes of rotation of the hinges that do form. 

FUTURE WORK 

GIRDS can be modified to handle a much broader base of 
problem than those which have been considered up to now. 
There are some modifications and additions that could be 
made to enhance the program’s performance. 

The most useful feature that could be added would be 
the capability to enter distributed loads. This task should 
not be too simplistic, however. Distributed plate loads do 
not transfer to the edges of the plate as linear line loads 
except in the case of circular plates and should not be 
handled in this way. Qne might approach this problem by 
linearizing or approximating the Fourier series solution for 
lcaded plates. 

Another area that improvement could be made is to allow 
for different yield moments in particular members. Adding 
arrays to keep track of this specific information and 
calling new variables as appropriate would be tedious but 


not particularly difficult from a programming standpoint. 
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APPENDIX A 


PROGRAM LISTING AND DIRECTIONS 





Grids is a very simple program to use. In the 
interactive mode it is self explanitory and requires no 
further information to run. To run the program more quickly 
when repetitive inputs are desired, an input data file may be 
used. Figure 19 shows the construction of such a file. One 
need only assign that file to be the source of input to run 
the program with all input read from the input file without 
further action from the user. 

The Grillage size limits placed on GRIDS was 
primarily for the convinence of the programmer. Arrays may 
be redimensioned to solve larger grids. The criteria 
provided in Appendix B for the IMSL linear programming 


routine must be met, however. 





SAMPLE INPUT 


=_— ss > 


6,1,592.2, 53. 3\ 
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Pee Soe. ©, 55.5 
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11,1,52.35,53.5 
di 65s)2,55.5 / 
ieteso cs. 55.35 / 
Wate ss. S, 55. 5/ 


o 


O. 


0,0, ¢ 


re 


There may be no blank 


blank for clarification. 


Lines in the data file. 


DATA FILE (WITH EXPANATIONS) 


GRID DIMENSIONS ON THIS LINE 
ROUNDARY CONDITIONS 
HORTIONTAL LENGTH 

VERTICAL LENGTH 

HORIZONTAL YIELD MOMENT 
VERFICAL YIELD MOMENT 


TYPE OF LOAD (1,2,0) 


REFERENCE NODE, DIRECTION (1.2), 
VALUE FOR POINT LOAD. OR A-NODE 
VALUE FOR LINE LOAD, LOCATION FOR 


POINT LOAD, OR R-NODE VALUE FOR LINE 


NO MORE LOADS OF THIS TFYPE FO ADD 
NO MORE LOADS TO ADD 


TYPE OF CONSTRAINTS (1,2,9) 


Lines were left 


There are four additional lines required 


if manual censtraints are used. 


Figure 19 Sample Input Data File 
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THIS BLOCK IS THE COMMON MEMORY SHARED BY ALL SUBROUTINES 
IT IS INCORPORATED BY USING THE INCLUDE ATATEMENT 


COMMON /A/ PTLOAD (20) ,NODTOT,NODESH, NODESV, 
II (30), Jd (30), VLEN, HLEN, XDCNST (20), 

LDIND (60), DIVAL (60) ,DILOC (40) ,LPT,LD2ND(60), 
D2AVL (60) , D2BVL (60) ,LNLDN, LDIS, NDCNST (20), 
LDZND (60) ,D2VAL (60) , IAUTO, NEQST 


+ + + + 


COMMON /B/ VMNT (S30) ,HMNT (20) , HORZMO, VERTMO 
COMMON /C/ ITOP, ILFT,IRGT, IRBTM, ITERAT 

COMMON /D/ PP, IER, IEQNM, NODEEG,NITRAT 

COMMON /E/ FOG(12,61),A2(41,61),B1 (12) ,B2(31) 
COMMON /F/ LIN(40) .MN, IA, LMAX,M1.M2, IW1, IRW1,LP2 


COMMON /GGG/ BHEMNT (20) , BHKLOC (20) , BVKEMNT (20) , BVKLOC (20) 


THIS PROGRAM SOLVES THE PROBLEM OF FINDING THE MAXIMUM LOAD FACTOR 
TO BE APPLIED TO A PARTICULAR LOADING CONDITION OF A RECTANGULAR GI 
IT SOLVES THE PROBLEM BY CONVERTING THE RESULTING SET OF EQUATIONS 
INTO A LINEAR PROGRAMMING PROBLEM. AT THE PRESENT TIME THE MAXIMU! 
SIZE OF GRID THAT CAN BE HANDLED IS A 6&6 X S GRID (IE 20 NODES ) 
THIS CONTRAINT I5 ARTIFICIAL AND IS DUE SOLELY TO THE VERY LARGE 
WORE VECTOR REQUIRED BY THE LP SUBROUTINE ZX3LP 

SEE IMSL LIBRARY REFERENCE MANUAL #8 FOR FURTHER INFOQ. 


DIMENSION A(82,29),B(85),C(29),PSOL (8%) .DSOL (82), 
+ RW(7200),I1W(240) 


FOR’ 


OPEN (UNIT=15, F ILE=* MOMENTS. DAT*® ,STATUS=’ NEW? ) 
PRINT#, “PROGRAM RESTRICTED TO A MAXIMUM GRID OF 6 X 3° 
PRINT#®,° ITE 20 NODES TOTAL OR LESS’ 
PRINT#, °° THIS RESTRICTION IS IMPOSED BECAUSE OF THE MAXIMUM’ 
PRINT*®.’SIZE OF THE WORK VECTOR USED BY ZXSLP IS CONSTRAINED’ 
PRINT*®, > BY THE PROGRAMMMER. ’ 
PRINT#, “SEE PROGRAM MANUAL FOR MORE INFOQ’ 


INPUT SIZE OF GRID TO BE SOLVED 
PRINT,’ 
PRINT*,” : 


PRINT*®.’ INPUT NUMBER OF NODES ALONG THE TOP AND SIDE OF THE GRID’ 
READ (S.*) NODESH, NODESV 

NODTOT=NODESH*NODESV 

IF (NODTOT.GT.3s0) GOTO 1 


INPUTCONDITIONS 








|. 


bh O 


A 


oon 


PRINT*,’ INPUT BOUNDARY CONDITIONS FOR THE TOP” 
PRINT*, 7LEFT,RIGHT AND BOTTOM EDGES OF THE GRID’ 
PRINT*,”’ IN THAT ORDER?’ 

PRINT*,’ INPUT A 1 FOR SIMPLY SUPORTED EDGE’ 
PRINT*, INPUT A 2 FOR A CLAMPED EDGE” 
PRINT#*, 7ALL VALUES ARE INTEGERS’ 
PRINT#, ’ ? 
PRINT*, ’ ? 
READ(S,*) ITOP,ILFT, IRGT, IBTM 
I=ITOP 
J=ILFT 
K=IRGT 
L=IBT 
IF (1.EQ@.1.0R.1.EQ.2.AND.J.EQ.1.0R.J.EQ.2. 
+AND.K.EQ.1.0R.K.EQ.2.AND.L.EQ.1.O0R.L. 
PRINT*, 7BOUNDARY CONDITIONS NOT ENTERED CORRECTLY. ° 
PRINT#,” TRY AGAIN’ 


GO TO 2 

CONTINUE 

PRINT*, 7 INPUT HORIZONTAL MEMBER LENGTH (REAL NUM. )’ 
PRINT#, ” ? 


READ(S,*) HLEN 
PRINT*,’ INPUT VERTICAL MEMBER LENGHT (REAL NUM.) ’ 
PRINT*,° , 
READ(S.*) VLEN 
PRINT*, °° INPUT HORIZONTAL BEAM YIELD MOMENT (REAL) ’ 
PRINT*, ’ ? 
READ(S5,*) HORZMO 
PRINT*,? INPUT VERTICAL BEAM YIELD MOMENT (REAL) ’ 
PRINT*, ° ? 
READ(5,*) VERTMO 
SET LOAD COUNTERS TO ZERO 
LPT=0 
LNLDN=0 
LDIS=0 
TEQNM=0 


NOW TO INITIALILZE THE LOAD AND MOMENTS ARRAYS 


CALL MINT 


CALL INPTL 


NOW TO ENTER ALL THE LOADS APPLIED TO THE GRID 


CALL ENTERL 


NOW TO ENTER ADDITIONAL CONSTRAINT INEQUALITIES 





7 © on 


i ial eat Ge 


es el alg 


15 


SEE USER’S HANDBOOK FOR FURTHER INFO 


PRINT*,’ THIS SUBROUTINE ALLOWS PLASTIC MOMENTS TO BE’ 
PRINT*, 7FORMED AT LOCATIONS OFF OF THE NODES OF THE’ 
PRINT*, ’GRILLAGE. THIS IS DONE BY WRITING ADDITIONAL’ 
PRINT*, ’CONSTRAINT EQUATIONS FOR THE LINEAR PROGRAMMNIG’ 
PRINT#*, 7 ROUTINE?’ 

PRINT*®,”7AUTO CONSTRAINTS WILL LOOK AT EACH MEMBER AND’ 
PRINT#, 7CONSTRAIN THE MAXIMUM MOMENT TO <= YIELD MOMENT’ 


PRINT#, ’ > 

PRINT, ” , 

PRINT*,* MANUAL ALLOWS THE USER TO CHOOSE WHICH MEMBERS’ 
PRINT*,”’ TO SET CONSTRAINTS FOR AND THEIR LOCATION ” 
PRINT#, ’ ? 

PRINT#, ’ ? 

PRINT*,”’ INPUT A 1 FOR AUTOMATIC CONSTRAINTS’ 

PRINT#, ” A 2 FOR MANUAL CONSTRAINTS’ 

PRINT#, ” A oO FOR NO CONSTRAINTS : 


READ(S,*) IAUTO 

IF (IAUTO.EQ.2) THEN 
CALL CONST 

ENDIF 


CONT INUE 


NOW TO CALCULATE THE NUMBER OF NODAL EQUATIONS 


NODEEQ= (NODESH—2) *# (NODESV—2) 


NOW TO CALCULATE THE NODAL EQUATIQN OF THE GRID AND WRITE THEM TO 
THE MATRIX FOG 


CALL EQNS 
WE MUST CALCUALTE THE DIMENSIONS OF ALL ARRAYS REFQRE 


FIXER BECAUSE ZXSLP REQUIRES IT 
NOW TO CALCULATE AND TRACK THE TOTAL NUMBER OF VARTIARLES 


IE TOTAL NUMBER OF NO-ZERO MQMENTS + LOAD FACTOR 


[=1 

DO 15, I1=1,NODTOT 
IF (HMNT(1).LT.0.5) GO TO 15 
San Ct) =T 
L=L+1 

CONTINUE 

LP2=L-1 

DO 16 I=1,NODTOT 
IF (VMNT(I).£T.0.5) GO TO 16 
LIN(L)=I 





L=L+1 


16 CONT INUE 


LMAX=L 


NOW TO MASSAGE THE INFORMATIGN CALCULATED AND STORED BY 
THE PROGRAM SO THAT IT CAN BE USED BY ZX3LP 


ele ee) 


Mi =LMAX+IEOQNM-1 
M2=NODEEQ 
IA=ML+M2+2 
IRWI=TAX TA+S¥EN1LA2EMN2+4 
IW1=2*M2+24M14+4 
MN=M1+M2 
C ADJUST M1 AND M2 IF THERE ARE EQUALITY CONSTRAINTS 
M1=M1—-NEQST 
M2=M2+NEQST 


CALL FIXER (A,B,C,PSOL, DSOL, RW, IW) 
C CALL BEAM CHECKING ROUTINE 
CALL BEAMCK 
ITERAT=ITERAT+1 
C CALL OUTPUT 
CALL OUTPUT 


BMAX=0.0 
DO 3799 JT=1,NODTOT 
IF (BHEMNT (1) .GT.HORZMO) THEN 
RMAX=BHEMNT (1) 
ENDIF 
IF (BVKMNT (1) .GT.VERTMO) THEN 
RMAX=BVKMNT (CT) 
ENDIF 
999 CONTINUE 
C LOOP TO WRITE EXTRA CONSTRAINT EQUATIONS IF AUTOMATIC IS SELECTED 
IF (TAUTO.EQ.1.AND. BMAX.GT.9.90. AND. 
+ ITERAT.LT.2) THEN 
CALL CONST 
GO TQ 5 
ENDIF 


END 


SUBROUTINE ENTERL 
INCLUDE * COMMON. FOR’ 


PRINT *.° LOADS ARE ENTERED REFERENCED TO THE NODE’ 








Brunt. 2. ABOVE OR TO THE LEFT OF THE LOAD’ 
PRINT +*, ALLOWED LOADING CONFIGURATIONS: ’ 

PRINT #, 1. POINT LOADS?’ 

PRINT +*,’ 2. LINEAR LINE LOADS ALONG A MEMBER’ 
Bilt ae, * =. LOAD ON A PLATE ELEMENT?’ 

Pint) ae, ” 

EmiMy *«,’ z 

PRINT *,”’ is 

Beit *,° ENTER THE FYFE @F LOAD TO ADD’ 

PRINT *, 7 1 = POINT LOADS’ 

PRINT *,?’ 2 = LINEAR LINE LOAD ALONG A MEMBER’ 
PRINT *,? QO = NO MORE LOADS TO ADD’ 


READ(S, *) MM 


IF (MM .EQ@. 9) GO TO 150 
IF (MM .EQ. 1) GO FO 100 
IF (MM .EQ. 2) GO TO 110 


CALL PNTLD (PTLOAD, NODTOT, NODESH, NODESV., VLEN 
+ ,HLEN,1II,JJ,LDIND, DIVAL, DiILOC,LPT) 


GO TO 90 

CALL LNLD (PTLOAD, NODTOT, NODESH, NODESV, VLEN 
+ ,HLEN,II,JJ,LD2ND, D2ZAVL, D2BVL,LNLDN) 

GO TO 90 

RETURN 

END 

SUBROUT INE PNTLD (PTLOAD.NODTOT, NQDESH, NODESV, VLEN 
+ »HLEN, II, JJ,LDIND,.DIVAL, DILOC.LPT) 


DIMENSIQN PTLOAD (NODTOT),DILOC(NODTOT*#2) 
DIMENSION DIVAL (NODTOT*2) 
DIMENSION II (NODTOT), JJ (NODTOT), LDIND (NODTOT*2) 


PRINT *, 7” THIS ROUTINE ADDS FOINT LOADS TO THE GRID. 
+ THERE MAY BE’ 

PRINT *.7A MAXIMUM OF ONE POINT LOAD PER NODE PLUS ONE 
+PER MEMBER?’ 

PRINT *, 7” THE REFERENCE NODE IS ABOVE OR TO THE LEFT OF 
+ THE LOAD’ 


PRINT *,° ENTER FOUR VALUES TO DEFINE EACH LOAD’ 
PRINT *,” 1ST VALUE = REFERENCE NODE (INTEGER)* 
PRINT *,7 ZND VALUE VERTICAL OR HORIZONTAL MEMBER* 





oe) 


PRINT *, 7’ 1 = HORIZONTAL MEMBER CINTEGER) ’ 
PRINT *,’ 2 = VERTICAL MEMBER CINTEGER) ’ 
PRINT *,’ QO = ON THE NODE CINTEGER) ’ 

PRINT *,’ SRD VALUE = LOAD MAGNITUDE (REAL) ’” 

PRINT *,’ 4TH VALUE = DISTANCE FROM THE REFERENCE’ 
PRINT *,”’ NODE (REAL) (IF THE 2ND VALUE = QQ’ 
PRINT *,”’ THEN THE 4TH MUST = QO)’ 


PRINT *,’ENTER ALL ZEROS WHEN COMPLETED WITH POINT LOADS’ 
READ (S,*) NDNUM, IDIR,VAL.RLOCA 


IF (NDNUM .E@. 9) GO TO 1530 
IF (NDNUM .GT. NODTOT) GO TO 190 
IF (II(NDNUM) .EQ. NODESV AND IDIR .EQ. 2) GOTO 1190 
IF (JJ CNDNUM) .EQ. NODESH AND IDIR .EQ. 1) GOTO 110 
IF (IDIR .E@. 2 .AND.RLOCA .GT. HLEN) GO TO 120 
IF (IDIR .EQ. 1 . AND. RLOCA .GT. VLEN) GO TO 120 
IF (IDIR .E@. 2) THEN 
RLOCA = -1.9 * RLOCA 
ENDIF 
LPT = LPT + 1 
MEGATIVE RLOCA INLICATES LOAD ON THE VERTICAL MEMBER 


LDIND(LPT) = NDNUM 
DIVAL(LPT) = VAL 
DILOC(LFT) = RLOCA 


ADD LOAD IF IT IS ON THE NODE 
IF (IDIR ~EQ. O ) THEN 
PTLOAD (NDNUM) = PTLOAD (NDNUM) +VAL 
GO TO 5 
ENDIF 


SET DUMMY VARIABLES FOR VERT OR HORIZ MEMBER AS APPROP 
IF (IDIR .~EQ@. 1) THEN 
NR=NDNUM+1 
RLEN = HLEN 
FLSE 
NB=NDNUM+NODESH 
RLEN = VLEN 
ENDIF 


CALCULATE AND ADD NODE LOADS FOR POINT LOAD APPLIED 
PTLOAD(NB) = PTLOAD(NB) + (CABS (RLOCA) *VAL/RLEN) 
PTLOAD(NDNUM) = PTLOAD  (NDNUM) + (RLEN-ARBS (RLOCA) ) *VAL/RLEN 


GO TO 3 


PRINT *,.° THE LAST NODE ENTERED IS NOT ON THE GRID,’ 
PRINT *, RE-ENTER ALL DATA FOR THE LAST LOAD’ 
GO TO 3 








PRINT *,”’ THE LAST LOAD ENTERED IS NOT ON THE GRID,° 
PRINT *, *RE-ENTER ALL DATA FOR THE LAST LOAD’ 
60 TO 


Of 


PRINT *,’LOAD LOCATION IS NOT ON THE MEMBER ADJACENT TQ’ 
PRINT *,” THE REFERENCE NODE, RE-ENTER ALL DATA FOR THE’ 
PRINT *,°LAST LOAD’ 

60 TO 3 


RETURN 
END 


SUBROUTINE LNLD(PTLOAD, NODTOT, NODESH, NODESV, VLEN 
+ ,HLEN, 1I,JJ,LD2ND.D2AVL, D2BVL,LNLDN) 


DIMENSION PTLCAD (NODTOT) ,D2ZAVL (NODTOT) , D2BVL (NODTOT) 
DIMENSION LD2NDCNODTOT),.IIC(NODTOT) , JJ (NODTOT) 


PRINT *.* THIS ROUTINE ADDS LINEAR LOADS. THE LOADS MUST BE 


+ON THE’ 
PRINT *,° GRID MEMBER. THE MAXIMUM NUMBER OF LINEAR LOADS 


+ ALLOWED’ 
PRINT *,”’ IS ONE PER GRID MEMBER. THE REFERENCE NODE IS 


+ AROVE OR’ 


PRINT #*,” TO THE LEFT OF THE LOAD. ‘ 
PRINT *,’ENTER FOUR VALUES TO DEFINE EACH LOAD ° 
PRINT *,’ IST VALUE = REFERENCE NODE (CINTEGER) 2 
PRINT *,’ 2ND VALUE = VERTICAL OR HORIZONTAL MEMBER i. 
PRINT *,’ 1 = HQRIZONTAL CINTEGER) : 
PRINT *,”’ 2 = VERTICAL CINTEGER) J 
PRINT #*,’ SRD VALUE = MAGNITUDE AT REFERENCE NODE ‘(REAL)’ 
PRINT *,’ 4TH VALUE = MAGNITUDE AT OPPOSITE NODE (REAL) * 


PRINT *, ENTER ZEROS FOR ALL FOUR VALUES WHEN COMPLETED 4 


READ (5,*) NDNUM, IDIR, AVAL, BVAL 


IF (NDNUM .E@. 90) GO TO 190 

IF (NDNUM .GT. NODTOT) GO TO 100 

IF (II (NDNUM) .EQ@. NODESV AND IDIR .EQ. 2) GOTO 1190 
IF (JJ¢NDNUM) .EQ@. NODESH AND IDIR .E@. 1) GOTO 110 
IF (IDIR .E@. 2) THEN 


RVAL = —-1.0*#BVAL 
AVAL=—1.Q#AVAL 
ENDIF 


LNLDN = LNLDN + 1 
NEGATIVE BVAL OR AVAL INDICATES LOAD ON THE VERTICAL MEMBER 








1190 


LD2ND(LNLDN) = NDNUM 
D2ZAVL (LNLDN) = AVAL 
D2BVL (LNLDN) = BYVAL 


SET DUMMY VARIABLES FOR VERT OR HORIZ MEMBER AS APPROP 
IF (IDIR .EQ. 1) THEN 


NR=NDNUM+1 
RLEN = HLEN 
ELSE 


NR=NDNUM+NODESH 
RLEN = VLEN 
ENDIF 
RVAL=ABS (BVAL ) 
AVAL =ABS5 (AVAL ) 
IF (AVAL .GT. BYVAL) THEN 


AA = 6.0 

BR = 3.90 

W = BVAL 

P = AVAL — BYVAL 
el OS) 

AA = 2.0 

BR = 6.0 

W= AVAL 

P = RBVAL —- AVAL 
ENDIF 


CALCULATE AND ADD NODE LCADS FOR LINE LOAD APPLIED 


PTLOAD (NDNUM) = PTLOADCNDNUM) +(¢ (WeHeRLEN/2) +P*RLEN/ BB) 
PTLOAD(NB) = PTLOAD(NR) + ((WeRLEN/2) + P*#RLEN/AA) 


60 TO 5 


PRINT *,* THE LAST NODE ENTERED IS NOT ON THE GRID,’ 
PRINT *, RE-ENTER ALL DATA FQR THE LAST LOAD’ 
GO TO 3 


PRINT *, °° THE LAST LOAD ENTERED IS NOT ON THE GRID,’ 
PRINT #,°* RE-ENTER ALL DATA FOR THE LAST LOAD’ 
GQ TO 5 


RETURN 
END 


SUBROUTINE MINT 


INCLUDE *~COMMON. FOR* 
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SUBROUTINE TO INITIALIZE THE BENDING MOMENT ARRAYS BASED ON THE 


GIVEN RBOUNDARY CONDITIONS 


K=0 

DO 11 I=1,NODESV 
DO 10 J=1,NODESH 
NDNUM=K*NODESH + J 
IF (ITOP.E@.1.AND.I.E£Q@.1) GO TO 100 
IF (ITOP.EQ@.2.AND.1I.£0.1) GO TO 150 
IF (IBTM.EQ.1.AND.I.EQ@.NODESV) GO TO 
IF (IBTM.EQ.2.AND.1I.EQ.NODESV) GO TO 
IF (ILFT.EQ.1.AND.J.EQ.1) GO TO 100 
IF (ILFT.EQ.2.AND.J.EQ@.1) GO TO =z00 
IF (IRGT.EG.1.AND.J.EQ.NODESH) GO TO 
IF (IRGT.EQ.2.AND.J.EQ.NODESH) GO TO 

VMNT (NDNUM) =1.0 

HMNT (NDNUM) =1.0 

TI (NDNUM) =I 

JJ (NDNUM) =J 


GO TO 10 
IF (J.E0.1.0R.J.E0.NODESH) THEN 
60 TO 100 
ELSE 
GO TO 200 
ENDIF 


VMNT C(NDNUM) =0.90 
HMNT (NDNUM) =0.0 
II (NDNUM) =I 
JJ (NDNUM) =J 
GO TO 19 
VMNT (NDNUM) =1.0 
HMNT (NDNUM) =0.0 
II (NDNUM) =I 
JJ (NDNUM) =J 
60 TO 190 
VMNT (NDNUM) =0.90 
HMNT (NDNUM) =1.0 
II (NDNUM) =I 
JJ (NDNUM) =J 
CONTINUE 
K=K+1 
CONTINUE 
RETURN 
END 


SUBROUTINE INPTL 


INCLUDE * COMMON. FOR’ 


100 
1590 


100 
500 
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THIS ROUTINE INITIALIZES ALL NODE LOADS TO ZERO 


DO 5S J=1,NODTOT 
PTLOAD (J) =0.0 
CONT INUE 

RETURN 

END 


THE FOLLOWING SUBROUTINE MASSAGES THE INFORMATION CALCUALTED AND 
STORED BY THE MAIN PROGRAM UNITS SO THAT IF CAN BE PASSED TO THE 
LINEAR PROGRAMMING SUBROUTINE ZX3LP 


SUBROUTINE FIXER (A,B.C.PSOL,DSOL,.RW, IW) 


DIMENSION ACTA.LMAX), BCIA) ,C(LMAX) ,PSOL (MN) ,.DSOL (IA), 
RWCIRWI). IWCTWI) 


INCLUDE ’*’ COMMON. FOR’ 
PART I--WRITE THE NODAL MOMENTS INEQUALITY CONSTRAINT EQNS 


DO 10 I=1,LMAX-1 
DO 20 J=1,LMAX 
IF(I-EQ@.J3) THEN 
AC(I.J)=1. 
ELSE 
A(I,J)=0. 
ENDIF 
CONTINUE 
IF(I.LE.LP2) THEN 
B(1)=HORZMO*2. 
ELSE 
R(1)=VERTMO*2. 
ENDIF 
CONT INUE 


NEXT WRITE THE ADDITIONAL CONSTRAINT EQUATIONS AS SELECTED BY THE 
USER IN THE LOAD ENTERING SEQUENCE 
IF (IEQNM.E@.9) GO TQ 100 
K=LMAXK+IEQNM-1 
DO 20 J=1,LMAX 

LL=1 

IF (J.E@.LMAX) THEN 
L=Z2*NODTOT+1 





OO) = oO 
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ELSE IF(J.LE-LP2) THEN 
L=LIN(J) 
ELSE 
L=LIN(J)+NODTOT 
ENDIF 
DO 40 I=LMAX,K 
ACI, J) =A2(LL,L) 
LL=LL+1 
CONTINUE 
CONTINUE 


LL=1 

DO SO I=LMAX.K 
B(I)=B2 (LL) 
LL=LL+1 

CONTINUE 

CONT INUE 


NOW TO WRITE THE NODAL EQUATIONS TO MATRIX A 


<=LMAX+IEQNM 
DO 40 J=1,LMAX 
Li=1 
IF (J.EQ.LMAX) THEN 
L=?*NODTOT+1 
ELSE IF(J.LE.LP2) THEN 
L=LIN(J) 
ELSE 
L=LIN(J)+NODTOT 
ENDIF 
DO 70 I=K,IA-2 
ACI,J)=FOG(LL,L) 
LL=LL+1 
CONTINUE 
CONTINUE 
LL=1 
DO SO I=K,IA-2 
B(I)=B1 (LL) 
LL=LL+1 
CONTINUE 


FINALLY WE WRITE THE C MATRIX 


DO 90 I=I,LMAX~1 
C(I)=0. 

CONT INUE 
C(LMAX)=1. 


NOW WE ARE READY TO CALL THE LP SUBROUTINE 
CALL ZXZLP (A, IA, B,C.LMAX,M1.M2,S5,PSOL, DSOL. RW, [W, TER) 








NOW TO MASSAGE OUTPUT OF ZXSLP SO THAT IT CAN BE USED 
IN THE MAIN PROGRAM OUTPUT ROUTINE 


es oo al la 


DO 99, I=1,LMAX-1 
K=LIN(I) 
IF(1.LE.LP2) THEN 
HMNT (K) =PSOL (1) -HORZMO 
ELSE 
VMNT (K) =PSOL (1) -VERTMO 
ENDIF 
99 CONTINUE 
PP=PSOL (LMAX) 


NOW TQ RETURN TO THE MAIN PROGRAM 


5 a Se al 


RE TURN 
END 


SUBROUTINE EQNS 
INCLUDE * COMMON. FOR’ 


SUBROUTINE TO CALCULATE AND STORE THE NODAL EQUATIONS OF THE GRID 


Get ©) 


FIRST STEP.—--INITIALIZE THE ARRAYS 
NODTOT=NODESH*NODESV 
DO 18 I=1,NODEEQ 
B1i(1)=0, 
DO 17 J=1,NODTOT*2+1 
FOG(1I,J)=0. 
17 CONTINUE 
18 CONT INUE 


C NCW TO CALCULATE THE NODAL EQUATIONS AND WRITE THEM TQ THE ROWS OF | 


IF (HLEN.GE.VLEN?) THEN 
XMUL T=HLEN 

ELSE 
XMULT=VLEN 

ENDIF 

M=NODESH+2 

N=NODTOT—-NQDESH-1 

L=1 


DO 19 I=M,N 
IF (JJ (1).EQ.NODESH.OR.JJ(1I).£Q.1) GO TO 19 
R=0. 
e CALCULATE THE COEFFICIENT OF HORZ MO. AT NODE I 
HC=2. /HLEN 
R=R+2. *HORZMO/HLEN 
FOG(L,1)=HC 
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CALCULATE THE COEFFICIENT OF HORZ MO. 


IF CINT(HMNT(I-1)).E@.0) GO TO 20 
HC=—1./HLEN 
R=R-HORZMO/HLEN 
FOG(L,1I-1)=HC 
CONTINUE 


TO 


CALCULATE THE COEFFICIENT OF THE HORZ MO. 


IF CINT CHMNT (I+1)).EQ@.9) GO TO 21 
HC=—-1./HLEN 
R=R-HORZMO/HLEN 
FOG(L,I+1)=HC 
CONT INUE 


CALCULATE THE COEFFICIENT OF THE VERT MO. 


VC=2. /VLEN 
R=R+2*VERTMO/VLEN 
K1=NODTOT+I 
FOG(L,K1)=VC 


CALCULATE THE COEFFICIENT OF THE VERT MO. 


K2=I-NODESH 
IF (INT (VMNT (K2)).EQ@.0) GO TO 22 
VC=-1./VLEN 
R=R-VERTMO/VLEN 
K3=NODTOT+K2 
FOG(L,K=)=VC 
CONTINUE 


CALCULATE THE COEFFICIENT OF THE VERT MO. 


K2=1+NODESH 
IF (INT (VMNT (K2)).EQ.9) GO TO 23 
VC=-1./VLEN 
R=R-VERTMO/VLEN 
KZ=NODTOT+K2 
FOG(L,K3)=VC 
CONT INUE 
K2=2*NODTOT+1 
FOG (L.K2)=-PTLOAD (I) 
B1(L)=R 


LEFT OF NODE I 


TO RIGHT OF NODE I 


AT NODE I 


ABOVE NODE I 


BELOW NODE I 


CHECK TO SEE IF RHS OF EQUATION IS NEGATIVE, BRANCH ACCORDINGLY 


IF(R.LT.O.) GO TO 24 
DO 25 J=1,2*NODTOT+1 

FOG (L,J)=FOG(L,J) *XMULT 
CONTINUE 

Ri (L)=B1 (L) *XMULT 
GO TO 26 
DO 27 J=1,NODTOT*2+1 

FOG (L,J)=FOG(L, J) *—-XMULT 
CONT INUE 

Bi (L) =B1 (L) *—XMULT 
CONTINUE 
L=L+1 
CONT INUE 
RETURN 





END 
SUBROUTINE CONST 


INCLUDE *’ COMMON. FOR’ 


THIS SUBROUTINE WRITES THE ADDITIONAL INEQUALITY CONSTRAINT 
EQUATIONS SELECTED BY THE USER TO ENSURE THAT THE BENDING 

MOMENT IN A HORIZONTAL OR VERTICAL BEAM DOES NOT EXCEED 

THE BEAM PLASTIC BENDING MOMENT (IE MOMENT BETWEEN NODES) 

THIS ALLOWS THE FORMATION OF PLASTIC HINGES AT LOCATIONS OTHER 
THAN THE NODES 

THIS ROUTINE ALSO HAS THE CAPABILITY TO CONSTRAIN MOMENTS OF A 
MEMBER AS HAVING ALREADY REACHED THE YIELD MOMENT. THAT IS 
SIMILAR TO HAVING A MEMBER THAT HAS ALREADY FAILED OR HAS MINIMAL 
STRENGTH LEFT DUE TO EXTREME LOCAL LOADING. THIS IS 

ACCOMPLISHED BY MAKING THE ADDITIONAL CONSTRAINT EQUATION 

AN EQUALITY CONSTRAINT RATHER THAN AN INEQUALITY CONSTRAINT 


GIG be) ba bala eae elects 


RE-INITIATE HMNT & VMNT ARRAYS FOR USE AS MULTIPLIERS TO 
BEQST=0 DETERMINE EXISTANCE OF PARTICULAR MOMENTS BASED ON BOUNDARIES 
5 CONTINUE 
10 IF (IAUTO.EQ.2) THEN 
PRINT*, ’ INPUT THE REFERENCE NODE FOR THE ADDITIONAL’ 
PRINT*, ’CONSTRAINT, AND A 1 FOR HORIZONTAL MEMBER’ 
PRINT*,°OR A 2 FOR VERTICAL MEMBER (INTEGERS) ’ 
PRINT*, ENTER ZEROS FOR NO CONSTRAINTS’ 
READ (5,¥*) NN,MANDIR 
IF (NN.EQ.0) GO TO 49a 
PRINT*, INPUT THE DISTANCE FROM THE REFERENCE NODE’ 
PRINT*,’TO LOCATE THE CONSTRAINT (REAL) ’” 
READ (5,%*) XX 
BHKLOC (NN) =XX 
BVKLOC (NN) =XX 
PRINT*, ENTER A 1 TO INDICATE THAT THE MEMBER HAS A’ 
PRINT*, "FULLY DEVELOPED YIELD MOMENT ALREADY DEVELCPED”’ 
PRINT*, NOTE: THIS TYPE MUST BE ENTERED LAST’ 
PRINT*, "ENTER A QO IF THE HINGE IS NOT ALREADY FORMED’ 
PRINT#*, ” (INTEGERS) ’ 
READ(5,*) N@ 
NEQST=NEQST+NQ 
GO TO 20 
ENDIF 


3 
ID 


DO 7900 NN=1, NODTOT 


BHPVAL=90.0 
BHPLOC=0.90 
RVPVAL=90.0Q 
BVPLOC=0.90 


J 
a) 


4 0 7 . = ‘ 








BHLAV=0.0 
BHLBV=0.0 
RVLAV=0.0 
BYLBV=0,0 


IF (TAUTO.EQ.2) THEN 
BHKMNT (NN) =HORZMO 
BYKMNT (NN) =VERTMO 


ENDIF 
DO 100 J=1,LPT 
IF (LDIND(J).EQ.NN.AND. DILOC(J) .£0.0.0) 
+ GO TO 100 
IF (LDIND(J).EQ.NN.AND. DILOC(J) .GT. 90.90 
= ~AND. BHKMNT (NN).GE.HORZMO) THEN 
RBHPVAL=DIVAL (J) 
RBHPLOC=D1LOC (J) 
ELSE IF (LDIND(J).EQ.NN.AND.DILOC(J).LT.Q.0 
a - AND. BVEMNT (NN) .GE.VERTMQ) THEN 
BVPVAL=ABS (DIVAL (J) ) 
BYPLOC=ABS (DILOC (J) ) 
ENDIF 
100 CONT INUE 
DO 200 J=1,LNLDN 
IF (LDOND(J). EQ.NN .AND. ((D2BVL(J) .GT. O.0O.O0OR. 
+ DZAVL(J).GT.9.9)) .~. AND. BHKMNT (NN) .GE.HORZMO) 
+ THEN 
BHLAV=D2AVL (J) 
RBHLBV=D2RBVL (J) 
ELSE IF (LD2ND(J).EQ.NN .AND. (D2BVL(J) .LT. Q. 
+ ~-OR. D2ZAVL (J).LT.90.9) . AND. BVKMNT (NN) 
+ -GE.VERTMO) THEN 
BVLAV=ABS (D2ZAVL (J) ) 
RVLBV=ARBS (D2RYVL (J) ) 
ENDIF 
200 CONTINUE 


C IF THERE ARE NO LOADS ASSOCIATED WITH THIS NODE, CONTNUE 


IF (BHPVAL.EQ.9.0 .AND. BYPVAL .EQ. 9.9 .AND. BHLAV .EQ. 











OM 


C 
C 


+ QO.0 .AND. BHLEBV .EQ. 0.0 .AND. BYVLAV .E@. 0.090 . AND. 
* RVLRBY .EQ. 9.9) GO TO 680 


IF (BHPVAL.EQG.0.AND. BHLAV.EQ.0.90. AND. BHLBV.EQG.0.0) 
+ GO TQ 450 
CHECK FOR VERTICAL OR HORIZONTAL FOR MANUAL MODE 
IF (MANDIR.EQ@.2) GO TO 450 
HORIZONTAL BEAM CALCULATIONS 
XL=HLEN 
A=RBHPLOC 
Q=BHPVAL 
STEP=HLEN/ 1000.0 
XX=BHKLOC (NN) 


INCREASING OR CONSTANT LINE LOAD WITH MAX MOMENT TO THE LEFT OF 
OR AT A POINT LOAD 
IF (BHLAV.LE.BHLBV.AND.XX.LE.A) THEN 
P=RHLRV—-BHLAV 
W=BHLAV 


X1=(1.-XX/XL) *HMNT (NN) 
2=(XX/XL) *HMNT (NN+1) 
2=(1.-A/ XL) #B#XX+. SH#WHEXLEXX+PHEXLEXX/G.0 


+ 2 DEWHX XBRL —-PHEXXHHS/ (6.O0%8XL) 
X4=HORZMO* (1.+(1.—-XX/XL) #HMNT CONN) + OXX/XL) 
+ *HMNT (NN+1) ) 
IDIR=1 


XDCNST (NN) =XX 

IF (NQ. EQ. 1) THEN 
NDCNST (NN) =—-NN 

ESE 
NDCNST (NN) =NN 

ENDIF 


CALL CONST2 (X1,X%2,.X5S,44, IDIR. NN) 


INCREASING OR CONSTANT LINE LOAD WITH MAX MOMENT 
TO RIGHT OF POINT LOAD OR NO POINT LOAD 
ELSE IF (BHLAV.LE.BHLBV.AND.XX.GT.A) THEN 


P=BHLBV-RBHLAV 
W=RHLAV 
XXX=XL—-XX 


X1=(XXX/ XL) *HMNT (NN) 

2=(1.-XXX/ XL) #HMNT ONN+1) 

X3=(A/XL) EQEXXX+. SHWHXLEXXX+PHXLEXXX/S. 
t—, SEWHEXXX¥#H#2 —. SHEPEXXXHH2ZtEPHXXXHHS/ (6.08XL) 








X4=HORZMO* (1.401. -XXX/XL) *HMNT ONN+1)4+¢XXX/XL) 


ts *HMNT CNN) ) 
IDIR=1 
XD IF (NQ.EQ.1) THEN 
NDENST (NN) =—NN 
ELSE 
NDCNST (NN) =NN 
ENDIF 


CALL CONST2 (X1,X%2,X%2,X4, IDIR, NN) 


C DECREASING LINE LOAD WITH MAX MOMENT TO THE LEFT OF OR AT 
C A POINT LOAD OR WITH NO POINT LOAD 


ELSE IF (BHLBV.LT. BHLAV. AND. XX.LE.A) THEN 
P=BHLAV—-BHLBV 
W=BRHLRBV 


X1=(1.—-XX/XL) #HMNT (NN) 
A2=(XX/ XL) *HMNT CNN+1) 
KS=C1.—-A/ XL) ROKK X+. DHWHXLEXX+PHXLEXK/ SS. — 2 GHEWEXKX HED 


+ +PHXXHHS/ (G6. O#XL) 
X4=HORZMO* (1.4401. —-XX/XL) *HMNT CONN) + OXX/XL) *HMNT (CNN+1) ) 
IDIR=1 


XDCNST (NN) =XX 
IF (NO. E®. 1) THEN 
NDCNST (NN) =—-NN 
ELSE 
NDCNST (NN) =NN 
ENDIF 


CALL CONST2 (X1,X2,X%3,X4, IDIR, NN) 


C DECREASING LINE LOAD WITH MAX MOMENT TO THE RIGHT OF THE 
C POINT LOAD 


ELSE IF (BHLBV.LT.BHLAV.AND.XX.GT.A) THEN 


P=BHLAV-BHVRV 
W=BHLRV 
XXX=XL—-X 


X1=¢€XXX/ XL) *HMNT CNN) 
X2=(1.—-XXX/XL) *HMNT CNN+1) 
X3=(A/KL) HBHXXX+. SHWHEXLEXXX+PHXLEXKK/ GS. —. DHWEKXX HHS 
+PHXXXH#HS/ (5. O8XL) 
X4=HORZMO*#(1.+¢01.—-XXX/XL) *HMNT (NN4+1) + ¢0XXX/XL) *HMNT CNN) ) 
IDIR=1 
XDCNST (NN) =XX 
IF (N@. EQ. 1) THEN 
NDCNST (NN) =—-NN 





EESE 
NDCNST (NN) =NN 
ENDIF 


CALL CONST2 (X1,X2,X2,X4, IDIR, NN) 
ENDIF 


C VERTICAL MEMBER CALCULATIONS 


430 CONTINUE 
C SKIP THIS SECTION IF THERE ARE NO LOADS ON THE VERT MEMBER 
IF (CBVPLOC.EQ.9.0.0R.BVPVAL.EQ.0.0).AND. BVLAV.EQ.0.0 
+ -AND. BVLBV.EQ@.9.9) GO TO S80 
C SKIP THIS SECTION IN MANUAL IF HORIZONTAL CONSTRAINT 
IF (MANDIR.E@.1) GO TO &80O 
XL=VLEN 
A=BVPLOC 
Q=BVPVAL 


XX=BVELOC (NN) 


C INCREASING LINE LQAD WITH MAX MOMENT TO THE LEFT OF 
C OR AT A POINT LOAD 
IF (BVLAV.LE.BVLBV.AND. XX.LE.A) THEN 
P=RVLBV—-BVLAV 
W=BVLAV 


X1=(1.—-XX/XL) *#VMNT CNN) 
a= (XX/XL) #VMNT ONN+NODESH) 
XS= (1. -—A/ XL) HQRXX+. DHWEXLEXX+PHXLEXX/ 56. —. DHWEX XHHL 


+ —P#XX##3/ (5. O#XL) 
X4=VERTMO* (1.4+¢1.-XX/XL) *VMNT (NN) + OXX/ XL) #VMNT ONN+NODESH) ) 
IDIR=2 


XDCNST (NN) =—-XX 
IF (NQ.EQ@.1) THEN 
NDCNST (NN) =—-NN 
EESE 
NDCNST (NN) =NN 
ENDIF 


CALL CONST2 (X1,X2,XZ,X4, IDIR, NN) 


C INCREASING LINE LOAD WITH MAX MOMENT TO RIGHT OF 
C POINT LOAD OR NO PT LOAD 
FLSE IF (BVLAV.LE. BVLBV. AND. XX.GT.A) THEN 


P=BVLBV—-RBVLAV 
W=BVLAV 
XXX=XL—-XX 











Xi=(XXX/XL) *#VMNT (NN) 
X2=(1.—-XXX/XL) *VMNT (NN+NOQDESH) 
KS=CA/SXL) KQHXXX+. SHWEXL HEX XX+PHXLEXXX/ Se — eo DHWHEXXX HED 


+ 2 DEP HEX XK HHZAEP HEX XX HHZ/ (CG. OHXL) 
X4=VERTMO* (1.+¢1.—-XXX/XL) ¥VUMNT (NN+NODESH) 

+ +(XXX/XL) *VMNT CNN) ) 
IDIR=2 


XDCNST (NN) =—-XX 

IF (N@Q.EQ.1) THEN 
NDCNST (NN) =—-NN 

ELSE 
NDCNST (NN) =NN 

ENDIF 


CALL CONST2 (X1,X2,X2,X4,IDIR.NN) 
C DECREASING LINE LOAD WITH MAX MOMENT TQ THE LEFT OF OR AT 
C A POINT LOAD OR WITH NO POINT LOAD 
ELSE IF (BVLBV.LT.BVLAV.AND.XX.LE.A) THEN 
P=RVLAV—-BVLREV 
W=RVL BV 


XL) #V¥MNT CNN+NOQDESH) 
S= C1. -A/ XL) FQEXX +. SEWEXLEXX+PHEXLEXK/ 3. — DEW EXX HED 


at; +tPHEXXH#HS/ (G.O8#XL) 
X4=VERTMO* (1. 4+¢01.-XX/XL) EVMNT CONN) + OXX/ XL) ¥VMNT CNN+NODESH) ) 
IDIR=2 


XDCNST (NN) =—XX 

IF (NQ.EQ.1) THEN 
NDCNST (NN) =—-NN 

ELSE 
NDCNST (NN) =NN 

ENDIF 


CALL CONST? (X1,X2,X2,X4, IDIR, NN) 


C DECREASING LINE LOAD WITH MAX MOMENT TO THE RIGHT OF THE 
C POINT LOAD 


ELSE IF (BVLBV.LT.BVLAV.AND.XX.GT.A) THEN 


P=BVLAV—-BVVBV 
W=RVLBV 
XXX=XL—-A 


X1=(XXX/XL) *VMNT CNN) 

X2=(1.—-XXX/XL) *VMNT (NN+NOQDESH) 

X= (A/SXL) HQEXXX+. DHEWHEXLEXXX*APHEXL EXXX/ Ge — 2 DEWEX XX HES 
+ +PHEXXXHHSS (G6.O#XL) 

X4=VERTMO*(1.+¢1.—-XXX/XL) #VMNT (NN+NODESH) 





Ag +OXXX/XL) *¥VMNT CNN) ) 
IDIR=2 

XDCNST (NN) =—-XX 

IF (N@Q. EQ. 1) THEN 
NDCNST (NN) =-NN 

BSE 
NDCNST (NN) =NN 

ENDIF 


CALL CONST2 (X1,X2,X3,X4, IDIR, NN) 


ENDIF 


580 IF (BHPVAL.EQ.0.0. AND. BVPVAL.EQ.90.90.AND. 


+ BHLAV.EFEQ.9.0. AND. BHLBV.ED.9O.0. AND. 
+ RBVLAV.EQ.0.0.AND. BVLBV.EQ.90.0.AND. 
+ TAUTO.EQ.2) THEN 
IF (MANDIR.EQ@. 1) THEN 
XMULT1=HMNT (NN) 
XMUL T2=HMNT (NN+1) 
XL=HLEN 
XMO=HORZMO 
IDIR=1 
ELSE IF (MANDIR.EQ. 2) THEN 
XMULTI=VMNT CNN) 
XMULT2=HMNT (NN+NODESH) 
XL=VLEN 
XMO=VERTMO 
IDIR=2 
ENDIF 
X1=¢€1.—-XX/XL) ¥®XMULTI1 
X2=(XX/ XL) *#XMULT2 
XS=0.0 


X4=XMO¥(1.+¢1.—-XX/XL) XMULTI+XX/XL*¥XMULT 2) 


IF (IDIR. EQ. 1) THEN 
XDCNST (NN) =XX 

ELSE IFCIDIR.EQ.2) THEN 
XDCNST (NN) =—-XX 

ENDIF 


IF (N@.EQ.1) THEN 
NDCNST (NN) =-NN 
ELSE 
NDCNST (NN) =NN 
ENDIF 
CALL CONST2(X1,X2.X3.X4, IDIR, NN) 
ENDIF 


690 IF (IAUTO.EQ.2) THEN 
PRINT*,’” , 


PRINT*, "ENTER 1 FOR ADDITIONAL CONSTRAINTS OR Q’ 
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PRINT#, FOR SQ ADDITICNAL RESTRAINTS REQUIRED’ 
READ (5, *) LI! 
IFC(IIIIT.£@.1) GO TQ 10 
GOTQ 720 
ENDIF 


CONTINUE 


RETURN 
END 
SUBROUTINE CONST2 (X1,X2,X%3,X4, IDIR, NDREF) 


INCLUDE ’ COMMON. FOR’ 
=2*NODTOT 
TEQNM=IEQNM+1 
2=NDREF +1 


IF(IDIR.EQ@.2) GO TO 20 
DO 10 I=1,NDEF-1 

AZ (IEQNM, I) =0. 
CONTINUE 

A? (IEQNM. NDREF) =X1 

A2 (IEQNM,N2) =X2 
DO 15 I=N2+1,L 

A? (IEQNM, I) =0. 
CONTINUE 

A2 (1TEQNM.L+1)=X3 

B2 (IEQNM) =Xx4 
GO TO 40 


CONTINUE 
J=NODTOT+NDREF 
<=.J+NODESH 
DO 25 I=1,J-1 

2 (IEQNM, I) =O. 
CONTINUE 

A2 (IEQNM,J)=X1 
DO ZO I=J+1,K-1 

AZ (IEQNM, I) =0. 
CONTINUE 

A2 (1EQNM.K) =X2 
DO 35 I=K+1,Lb 

A2(IEQGNM,1)=0. 
CONTINUE 

AZ (IEQNM.L+1) =X 

B2 (IEQNM) =x4 
CONTINUE 


RETURN 
END 








AMES=" BHLBYV’ 


999 


2900 


SUBROUTINE BEAMCK 


INCLUDE ’ COMMON. FOR’ 
CHARACTER* 9 NAME. NAME1, NAME2, NAME2, NAME4, NAMES 
CHARACTER® 9 NAMES, NAME7, NAMES 


NAME 1=’ I/NODE#’ 
NAME="LD2ND #° 
NAME2=’ DZAVL’ 
NAME2=’ D2RB 


NAMES=" BYLAYV’ 
NAME7=’ BVLBV’ 
NAMES=’ NNN/LNLDN* 


DO 999 [T=1,LNLDN 
CONTINUE 


DO 700 I=1, NODTOT 


RBHEMNT (1)=0.0 
BHKLOC (1)=0.90 
RVEMNT (1) =0.0 
RBYKLOC (1) =0.0 
BHPVAL=0.0 
BHPLOC=0.0 
BYVPVAL=0.0 
RVPLOC=0.90 
BHLAV=0.90 
BHLRBYV=0.0 
RVLAV=0.90 
RVLBV=0.Q0 


FORMAT (2X, A9,1X,14,35X, au: iat ll AY, G12.2, 5X, AP, 


DO 100 J=1,LPT 


IF (LDIND(J).EQ.I1.AND. DILOC(J) .£@.0.0) GO TO 100 
IF (LDIND(J>).EQ.1I.AND. DILOC(J) .GT. 9.0) THEN 


BHPVAL=D1VAL (J) #PP 
BHPLOC=D1iLOC (J) 


FLSE IF (LDIND(J).EQ@.1.AND.DILOC(J).LT.9.9) 





RVPVAL=ABS (DI1VAL (J) ) #PP 
BVPLOC=ABS (DILOC (J) ) 


ENDIF 


100 CONT INUE 


DO 200 NNN=1,LNLDN 
IF ((LD2ND(NNN).EQ@.1) .AND. (D2BVL (NNN) .GT. 0.0.0R. 
& D2AVL (NNN) .GT.0.0)) THEN 


RHLAV=D2AVL (NNN) *#PP 
RBHLBV=D2BVL (NNN) #PP 


ELSE IFC (LD2ND(NNN).EQ.I1). AND. (D2BVL (NNN) .LT.0.90.O0R. 
+ D2ZAVL (NNN) .LT.9.9)) THEN 


RVLAV=ABS (DZAVL (NNN) ) ¥PP 
RVLBV=ABS (D2BVL (NNN) ) PP 


ENDIF 


hJ 


Q0 CONT INUE 
C IF THERE ARE NO LOADS ASSOCIATED WITH THIS NODE, CONTNUE 


IF (BHPVAL.E@Q.0.0 .AND. BYVPVAL .EQ. 9.90 .AND. BHLAV .EQ. 
ae 0.0 .AND. BHLBV .FQ. 90.0 .AND. BVLAV .EQ. 9.0 .AND. 
+ BVLBV .E@. 9.9) 60 TO 450 


IF (BHPVAL.EQ.0. AND. BHLAV.EQ.0.0. AND. BHLBV.EQ.9.Q0) 
+ 60 TO 4590 


C HORIZONTAL BEAM CALCULATIONS 
XL=HLEN 
=BHPLOC 
Q=BHPVAL 
BM1=HMNT (I) 
BM2=HMNT (I +1) 
STEP=HLEN/19000.0 


BMXA=0.0 
RBMXB=0.0 
BMX=90.0 


DQ 300 NNN=1, 1000 
X=NNN*STEP 





XX=X 
C SET X FOR SECOND INTERVAL IF TO RIGHT QF POINT LOAD 
IF (A.NE.O.0.AND.Q.NE.0.90.AND.X.GT.A) THEN 
X=X—-A 
END IF 


C INCREASING OR CONSTANT LINE LOAD 
IF (BHLAV.LE.BHLRBV) THEN 
P=BHLBV—-RBRHLAV 


W=RHLAV 
RCOM= (BM2—-BM1) #X/XL4+WHXL*X*. 5+ 
+ P*XL#X/6.0-. DHWHEX HHL —PHEX HRS / 
2 (6.0%XL) 
C DECREASING LINE LOAD 
ELSE 
P=BHLAV—-BHLBV 
W=RHLEV 
RCOM= (BM2—-RM1) *#X/XLAWHEXLEX#. 5+ 
at P#XL#EX/ 2. O-. DEWEX HHL -—PHEXHHLELS 
+ +P#XHHS/ (5. O#XL) 
ENDIF 


IF (Q.£6Q.0.0.0R.A.E8@.9.90) THEN 
RMX=RCOM+BM1 

ELSE IF (X.LE.A) THEN 
RMXA=RCOM+BM1+ (1-A/ XL) #O*#X 


ELSE 
RMX B=BMXA+RCOM—-A#Q#X / XL 
ENDIF 
IF (ABS (RMX) .LT.~ABS (BMXA).QR.ABS (BMX) 
+ -LT.ABS(BMXB)) THEN 
IF (ABS (BMXA).GT.ABS(RBMXR)) THEN 
RMX=BMXA 
ELSE 
RMX=RMXR 
ENDIF 
ENDIF 


IF (ABS (BHKMNT(I)).LT.ABS(BMX)) THEN 
BHKMNT (1) =BMX 
BHKLOC (1) =XxX 


ENDIF 
300 CONTINUE 
4590 CONTINUE 
C SET MAX MOMENT TO CONSTRAINED VALUE IF MEMBER HAS EXTRA 
C CONSTRAINT EQUATION FOR IT AND CONSTRAINT 15 LIMITING 


IF (XDCNST (1).GT.90.90. AND. (BHKMNT (I) .GE.HORZMO 





- »«OR.NDCNST (I).LT.0) ) THEN 
RHEMNT (I) =HORZMO 
BHKLOC (1) =XDCNST (I) 


ENDIF 
C VERTICAL MEMBER CALCULATIONS 
C SKIP THIS SECTION IF THERE ARE NO LOADS ON THE VERT MEMBER 


IF ( (BVPLOC.E@Q.0.0R. BVPVAL.EQ.0.90).AND. (BVLAV.EQ.0.0 
- AND. RVLRV.E@.9.0)) GO TO 610 


A=BVPLOC 

Q=BVPVAL 
BMI=VMNT (7) 
RM2=VMNT CI+NODESH) 
STEP=VLEN/1000.0 
XL=VLEN 


BMXA=0.0 
BMXR=0.0 
EMX=0.0 


DO 600 NNN=1, 1900, 
¥=NNN#*STEP 
XX=X 

C SET X FOR SECOND INTERVAL IF TO RIGHT OF POINT LOAD 
IF (A.NE.0.0. AND. Q.NE.0.0.AND.X.GT.A) THEN 


Cc INCREASING OR CONSTANT LINE LOAD 
IF (RVLAV.LE.BVLBV) THEN 
F=RVLBV-BVLAV 
W=BRVLAV 


RCOM= (BM2—-BM1) © (X/ XL) tWHEXLEX*®. O+ 


+ (P¥#XL#¥X/6.0) —(. SEWHK HH) — 
+ (P¥X##5/ (5. O#XL) ) 
C DECREASING LINE LOAD 
ELSE 
P=BVLAV—-RVLRBV 
W=RVLRV 
RCOM= (BM2—-BM1) #X/XL+WHXLEX*. OF 
+ P*XL#X/3.0—, SHEWHEXHHL—PHXHHL ELS 
+ +tPEXHHS/ (5. O#XL) 


ENDIF 


IF (Q.EQ.0.0.0R.A.EQ.90.0) THEN 
BMX=RCOM+BM1 
ELSE IF (X.LE.A) THEN 


>» 
> 


a 





BMXA=RCOM+BM1+ (1-A/XL) *O*X 


SESE 
RMXB=BMXA+RCOM—-A*#O*X / XL 
ENDIF 
IF (ABS (BMX) .LT.ABS (BMXA) .OR. ABS (BMX) 
a »~-LT.ABS (BMXB)) THEN 
IF (ABS (BMXA).GT.ABS (BMXB)) THEN 
RBMX=RMXA 
ELSE 
BMX=BMXB 
ENDIF 
ENDIF 


IF (ABS (BVEMNT (I) ).LT.ABS(BMX)) THEN 
BVKMNT (I) =BMX 
RBVKLOC (I) =X 


ENDIF 
600 CONTINUE 
610 CONTINUE 


IF (XDCONST(I).LT.9.0 .AND. (BYVKMNTCI) .GE. VERTMO 
at »~OR. NDCNSTC(I).LT.9.0)) THEN 
RVEKMNT (I) =VERTMO 
BVKLOC (1) =ABS (XDCNST (T) ) 
ENDIF 


790 CONTINUE 


RETURN 
END 


SUBROUTINE OUTPUT 
INCLUDE *COMMOQN. FOR’ 


CHARACTER*#15 SIDE. BBCC, DIR, FLAG 
CHARACTER*¥20 CNST 


C THIS ROUTINE OUTPUTS THE GRID LOADING SCHEME AND ALL MOMENTS AND LO. 
C FACTORS 
C CHECK ITERATION NUMBER IF IT IS ’17 THEN PRINT GRID DATA 
IF CITERAT.NE.1) GO TO SO 
C PRINT OUT GRID SIZE 


WRITE (15,1000) 
WRITE (15,1100) NODESH.HLEN 
WRITE (15,1200) NODESV, VLEN 


C PRINT QUT THE GRID BOUNDARY CONDITIONS 





WRITE (15,1200) 
DQ HEN 
IBOUND=ITOP 
SIDE=’ TOP’ 
FLSE IF (1.EQ@.2) THEN 
IBOUND=IBTM 
SIDE=’ BOTTOM’ 
FLSE IF (1.EQ.3) THEN 
ITROUND=ILFT 
OIDE="LEFT’ 
FLSE 
IBOUND=IRGT 
SIDE=? RIGHT’ 
ENDIF 


IF (ITBOUND .EQ. 1) THEN 
BPCC=’SIMPLE SUPORTED’ 
FLSE 
RBCC=’* CLAMPED’ 
ENDIF 


WRITE (15,1400) SIDE. BBCC 


15 CONTINUE 
PRINT*#, HORZMO,. VERTMO 

C WRITE THE ZX=LP ERROR NUMBER 
WRITE (15,1425) HORZMNO, VERTMO 


1425 FORMAT (////,” HORZ BEAM PLASTIC BENDING MOMENT = °*,G612.2, 
+//,”* VERT BEAM PLASTIC BENDING MOMENT = *,612.2 


WRITE (15, 1450) TIER 


1459 FORMAT (/7/,5X,” ZX2LP ERROR NUMBER = *,14) 
C PRINT OUT ALL THE POINT LOADS ON THE GRID 
WRITE (15,1500) 
WRITE (15, 1500) 


DO 25 I=1,LPT 
IF (ABS (DI1LOC(1I)).LT.0.01) THEN 
DIR=’ON THE NODE’ 


ELSE IF (D1ILOC DIR=*" TQ THE RIGHT’ 
ELSE 

DIR="* BELOW?’ 
ENDIF 


WRITES (15,1700) LDIND(I), DIR, ABS(DILOC(I)), DIVAL (I), 1 
20 CONTINUE 


ie PRINT QUT LINE LOADS 
WRITE (15,1800) 











DO 30 I=1,LNLDN 
IF (D2RVL(I) .GT. 0.0) THEN 
NB=LD2ND (I) +1 
ELSE 
NB=LD2ND (1) +NQDESH 
ENDIF 
WRITE (15,1900) LD2ND(I) NB, ABS(D2ZAVL(I)), ABS(D2ZBVL(1I)), 1 
CONTINUE 


FORMAT (’17,5X, GRID CHARACTERISTICS AND LOADING’ ,///) 
FORMAT (6X, NUMBER OF NODES HORIZONTAL = ’,13,5X, 
+7HORIZONTAL LENGTH =’ ,F4.2) 
FORMAT (&X,’*NUMBER OF NODES VERTICAL = ’,13,5%, 
+? VERTICAL LENGTH =’,F4.2) 
FORMAT (’0’,5X, BOUNDARY CONDITIONS’. /) 
FORMAT (6X,06,725,A15) 
FORMAT (///,6X,”POINT LOADS ADDED TO THE GRID’,/) 
FORMAT (3X, REF NODE’,T14,’DIR FM NODE’ ,T25, 
+’ DISTANCE FROM NODE’, 157, MAGNITUDE’ ,T79, 
+7POINT LOAD NUMBER? ) 
FORMAT (6X,1I3, 714,013,745, FS. 2, 158, F 10.2, 172, 13) 
FORMAT (8¢/),5X,’LINE LOADS’,//,6X,” BETWEEN NODES’, 
+T=0, REF NODE MAG’,T45, OTHER NODE MAG’, 
+T65, 7LINE LOAD NUMBER’ ,/) 
FORMAT (6X,I3,”’ AND *, I3,728,F10.2,1T44,F10.2,770, 15) 


CONT INUE 


PRINT QUT THE FINAL NODAL LOADING THAT CALCULATIONS ARE BASED ON 


WRITE (15, 2200) 
FORMAT (8¢/7),5X,* NODE LOADS USED IN CALCULATIONS’ ,//) 
WRITE (15, 2400) 
FORMAT (6X, 7NQDE NUMBER’ ,T20, ’NQDE LOAD’) 
DO SO I[=1,NODTOT 


WRITE (15,2500) I,PTLOAD(I) 
FORMAT (10X, 13, T21,F8.2) 


CONT INUE 
WRITE (15, 2550) ITERAT 


PRINT OUT THE EXTRA CONSTRAINT DATA 
WRITE (15,2525) 
FORMAT (//,2X, "ADDITIONAL CONSTRAINTS USED’,/) 


IF (IAUTO.EQ.1) THEN 
WRITE (15, 2527) 
ELSE IF (IAUTO.EQ. 2) THEN 





Fae ae A 
2u20 
a 7 


= 


WRITE (15, 2528) 
ELSE 

WRITE (15,2529) 
ENDIF 
FORMAT (2X, 7 AUTOMATIC CONSTRAINTS WERE USED’ ) 
FORMAT (2X, 7MANUAL CONSTRAINTS WERE USED’) 
FORMAT (2X. ’NO ADDITIONAL CONSTRAINTS WERE USED’ ) 


WRITE (C15, 2540) 
FORMAT (6X, ’BETWEEN NODES’ ,2xX, ’LOCATION’, 


3a : TYPE USED’ ) 


DO 57, I=1,NODTOT 
NOD1=ABS (NDCNST (I) ) 

IF (XDCNST(1).GT.9.0) THEN 
NOD2=ABS (NDCNST(I)) +1 

ELSE IF (XDCNST(1I).LT.9.9) THEN 
NOD2=ABS (NDCNST (I) ) +NODESH 

ELSE 

GO TO 57 

ENDIF 

IF (NDCNST(1I).GT.0.0) THEN 
CNST=" INEQUALITY ” 

ELSE 
CNST=" EQUALITY ” 

ENDIF 


WRITE (15, 2545) NOD1, NOD2, ABS (XDCNST (1) ) ,CNST 
CONTINUE 
FORMAT (6X,1IZ,” AND °,12,6X.G10. 2, A20) 

PRINT OUT MOMENTS AT EACH NODE AND THE LOAD FACTOR 


FORMAT (2%,//, 7 *#*** ITERATION NUMBER IS’,14,° “*#*#’.//) 
WRITE (15,2400) 
FORMAT (2(/).6X,’MOMENTS AND LOAD FACTOR’) 
WRITE (15,2790) 
FORMAT (//,6X,’ NODE NUMBER’ ,T25, "HORIZ MNT’,1T40,° VERT MNT”) 


DO 60 I=1,NODTOT 
WRITE (15, 2800) I, HMNT (I), VMNT (I) 

FORMAT (8X, 13,1726,ES8. 2, T42,£8. 2) 

CONTINUE 

WRITE (15.2900) PP 

FORMAT (///,6X,* THE GRID LOAD FACTOR = ’,G8.2) 


PRINT QUT MAX MOMENT AND LOCATION FOR EACH MEMBER 





WRITE (15, 3000) 
3000 FORMAT (3(/),////,6X%. MAXIMUM MOMENTS AND LOCATION 
+ ON LOADED MEMBERS’, /) 
WRITE (15,4050) 
4050 FORMAT (6X, 7 BETWEEN NODES’, 130, ’MAGNITUDE’ , 145, 
+ ’*LOCATION?) 


WRITE (15,4100) 
4100 FORMAT (6X, "HORIZONTAL MEMBERS?’ ) 
DQ 70 I=1,NODTOT 
IF (BHEKMNT(1I).EQ.0.0) GO TO 70 
IF (ABS (BHKMNT (I)).GT. (HORZMO+.01)) THEN 
FLAG=’ #*? 
ELSE 
FLAG=" ’ 
ENDIF 
WRITE (15,4200) I,1+1,BHKMNT(I),BHKLOC(I),FLAG 
79 CONTINUE 


WRITE (15,4300) 
4300 FORMAT (6X, ’ VERTICAL MEMBERS’ ) 
DO 80 I=1,NODTOT 
IF (BVKMNT (I).EQ.0.0) GO TQ gO 
IF (ABS (BVKMNT(1I)).GT. (VERTMO+.01)) THEN 
FLAG =’ «*’ 
ELSE 
FLAG =" ” 
ENDIF 
WRITE (15,4200) 1, I+NODESH, BVKMNT (I) , BVKLOC (I) , FLAG 
B80 CONTINUE 
WRITE (15,5000) 
FORMAT (6X,1I3,’ AND °,1%,728.G610.4,8X,G610.4,2X,A15) 
FORMAT (////,5%,’** INDICATES > YIELD MOMENT") 
RETURN 
END 


ie i 


LA > 
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APPENDIX B 


IMSL LINEAR PROGRAMMING ROUTINE 








The linear programming routine used by GRIDS is taken 
from the International Mathematics and Statistics Library (IMSL). 
GRIDS uses routine ZX3LP. The attached pages describe the 
routine input requirements and output formatting and variables. 
This information is pravided for a better understanding of the 
dimensioning process that GRIDS performs. These pages have been 


reproduced from the IMSL library reference manual. 





PURPOSE 


USAGE 


ARGUMENTS 


F 


TA 


DSOL —- 


RW 


| 


Solve the linear programming problem via the revise 
simplex algorithm —- easy to use version 


CALL ZXZLP (A,IA,B,C,N,M1,M2,S,PSOL, DSOL, RW. IW, IER) 


Matrix of dimension M14+M2+2 by N containing the 
coefficients ot the Ml rows followed by the 
coefficients of the M2 equality constraints (input). 
The last two rows are used only as working storage. 


Row dimension of matrix A exactly as specified 
in the dimension statement in the calling 
program. (input) Two rows of A are required for 
working storage, and therefore, [A must not be 
less than M1+M2+2. 


Vector of length M1+M2+2 containing the right 
hand sides of the inequality constraints. (input) 
The last two elements of B are used as working 
storage. 


Vector of length N containing the coefficients 
of the objective function. (input) 


Number af unknowns in the model. (input) 
Number of inequality constraints. (input) 
Number of equality constraints. (input) 
Value of the objective function. (output) 


Vector of length N containing the primal 
solution. (output) PSOL is also used as work 
storage and therefore must have length at least 
MAX (n, M1+M2). 


Vector of length M1I+M2+2 containing the DUAL 
solution. (output) That is, DSOL (1), ..., 

DSOL (M1+M2) contain the solution to the problem 
MIN BT*¥Y subject to AT*Y is greater than or 
equal to C and Y greater than or equal to OQ 
where AT = A-transpose and BT = B-transpose. 
When the primal problem has equality 
constraints, the corresponding components of the 
dual solution are unconstrained. DSOL (M1+M2+1> 
and DSOL (M1+M2+2) are used as working storage. 


Work Vector of length (M14+M2+2)*(ML+M2+2) + S*M1 
Eee eats. 
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IW - Work Vector of length 2*M2 + 3*M1 + 4. 


IER - ERROR indicator. (output) 
Terminal Error 
IER = 120 Indicates that IA is less than 
Mirher2. 
IER = 131 Indicates that the cost criterion 
has unbounded values. 


IER = 122 Indicated that the Maximum number 
of iterations was reached in ZXOLP. 
IER = 133 Indicated that no feasible solution 


exists. 
Warniq (with fix) 

IER = 79 Indicates that some artifical 
variables remain in the solution 
basis at a zero level after phase 
i This condition can be caused by 
having redundant constraints. 
Nevertheless, a solution is computed 
and returned in PSOL and DSOL. 


Algorithm 


To solve the linear programming problem, 


Maximize C,PsoL, se ree + CyPSoL,, = SP 
subject to 
a, ,PsoL, fo tenet a. PSOL,, < RB. 1=1,..-5M, 
a ,PsoL, ee ene = a. yPsol,, = BR. tr oe 
Foie Zoe jrigecs 
where M = Mil + M2. 
Fhe DUAL linear programming problem is, 
Minimize b,DSOL, ea ca + BYyVSol,, = $D 
subject to 
a ess - Py ee 
Stee ai oye eler = a j : 
DSOL . > O T=)... 6 = et 
1 sa 1 


DSOL . unrtestricted in sign when i = Mi tise ee aM 





ZXSLP computes the solution to the primal problem, PSOL. the 
solution to the dual problem, DSOL. and the values of the 
objective function S=SP=SD. 


ZXSLP calls ZXOLP which solve the linear programming problem 
by the revised simplex method. 


where M = M1 + M2. 
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APPENDIX C 


SAMPLE PROBLEMS AND RESULTS 


=o) = 


EXAMPLE 1 


GRID CHARACTERISTICS AND LOADING 


NUMBER OF NQDES HORIZONTAL = 6 HORIZONTAL LENGTH = 10.900 
NUMBER OF NODES VERTICAL = a VERTICAL LENGTH = 19.00 


ROUNDARY CONDITIONS 


TOP CLAMPED 
ROTTOM CLAMPED 
LEFT CLAMPED 
RIGHT CLAMPED 


HORZ BEAM PLASTIC BENDING MOMENT = 0, LOE+O5 
VERT BEAM PLASTIC BENDING MOMENT = 0, 1OE+05 
ZX35LF ERROR NUMBER = QO 
POINT LOADS ADDED TQ THE GRID 
REF NODE DIR FM NODE DISTANCE FROM NODE MAGNITUDE LOAD NUMBER 
8 ON THE NODE 0,00 1.00 1 


LINE LOADS 


RETWEEN NODES REF NODE MAG QTHER NODE MAG LINE LOAD NUMBER 
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NODE LOADS USED IN CALCULATIONS 


NODE NUMBER NODE LOAD 
1 Q.00 
2 0.00 
Zz 0.00 
4 0.00 
5 0.00 
& 0.00 
7 0.00 
8 1.00 
9 Q.00 

10 0.00 
at 0.00 
Fe 0.00 
1S Q.00 
14 2.00 
15 0.00 
16 0.00 
17 0.00 
18 0.00 
19 0.00 
20 0.00 
21 0.00 
22 0.00 
DS 0.00 
2 0.00 
25 0.00 
26 0.00 
2 Q.00 
2 9.00 
2 0.00 
=a Q.00 


##*#* ITERATION NUMBER IS 1 KEE 


ADDITIGNAL CONSTRAINTS USED 


NO ADDITIONAL CONSTRAINTS WERE USED 
RETWEEN NODES LOCATION TYPE USED 


MOMENTS AND LOAD FACTOR 





= 


NODE NUMBER HORIZ MNT VERT MNT 
1 0. 00E+00 0. O0E+00 
2 0. 00E+00 —. 10E+05 
a 0.G0E+00 —. LOE+OS 
4 0.90E+00 —~. LOE +06 
oa Q. 00E+00 —. 10E+05 
6 0. 00E+00 0. O0E+00 
7 =. 10E+¢3 0. O00E+00 
8 0. 10E+05 0.10E+05 
9? —.10E+05 0. 20E+04 
—. 40E+04 —. 1OE+05 
—. 70E+04 —. 10E+05 
—. 10E+05 0. 00E+00 
—. 1Q0E+05 0. Q0E+900 
0. 84E+04 =. 102105 
0. 68E+04 —.10E&+05 
—. 1Q0E+05 —~. 10E+04 
—.89E+04 —.10E+05 
=. 1OE+05 0. 00E+00 
—.10E+05 0. 00E+00 
—.10E+05 —« 1623-05 
—. 10E+05 —. 7/3E+04 
—.O0E+04 =. 1OE2O> 
—. 19E+05 —. /8E+04 
ae OE +OS 0. 00E+00 
0. OOE+00 0. 0OOE+00 
Q. O0E+00 —, 1OE+05 
0.Q0E+00 —. 10E+05 
0. Q0E+00 —. 1OE+05 
0, OOE+00 —. 1O0E+05 
0. 00E+00 0. O0E+00 


THE GRID LOAD FACTOR = 0.80E+04 


MAXIMUM MOMENTS AND LOCATION ON LOADED MEMBERS 


RETWEEN NODES MAGNITUDE LOCATION 
HORIZONTAL MEMBERS 


VERTICAL MEMBERS 


*%* INDICATES > YIELD MOMENT 
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EXAMPLE 2 


GRID CHARACTERISTICS AND LOADING 


NUMBER OF NODES HORIZONTAL = oO HORIZONTAL LENGTH = 190.900 
NUMBER OF NODES VERTICAL = =) VERTICAL LENGTH = 10.900 


RQUNDARY CONDITIONS 


TOP SIMPLE SUPORTED 
BOTTOM SIMPLE SUPORTED 
Per SIMPLE SUPORTED 
RIGHT SIMPLE SUPORTED 


HORZ BEAM PLASTIC BENDING MOMENT = 0O.10E+05 
VERT BEAM PLASTIC BENDING MQMENT = 0. 10E+05 
ZX=LP ERROR NUMBER = QO 
POINT LOADS ADDED TO THE GRID 
REF NODE DIR FM NODE DISTANCE FROM NODE MAGNITUDE LOAD NUMBER 
LINE LQADS 
BETWEEN NODES REF NODE MAG OTHER NODE MAG LINE LOAD NUMBER 
7 AND o 1.00 1.00 1 
7 AND Pia 1.00 1.00 = 
8 AND 12 1.00 Hele. = 
12 AND > 1.00 1.00 2 





a a= 


NODE LOADS USED IN CALCULATIONS 


NODE NUMBER NODE LOAD 


1 0,00 
ae 0.00 
3 0,00 
4 0.00 
a 0,00 
& 0.00 
7 10.00 
8 10.00 
9 0.00 
10 0.00 
11 0,00 
ee 10.00 
13 10, 00 
14 0.00 
15 0.00 
16 0.00 
a7 0.00 
18 0.00 
7 0,00 
20 0,00 
z= 0.00 
ee O.Q0 
Fe 0.00 
Pia 0.00 
au 0.00 


*#*e%% ITERATION NUMBER IS 1 HEH 


ADDITIONAL CONSTRAINTS USED 


AUTOMATIC CONSTRAINTS WERE USED 


RE TWEEN NODES LOCATION TYPE USED 

MOMENTS AND LOAD FACTOR 

NODE NUMBER HORIZ MNT VERT MNT 
1 0, OOE+00 QO, OOE+00 
= 0, OOE+00 0. QODE+00 
> 0, 00E+00 0. 90E+00 
a O, OOE+00 0, QOE+00 
oa 0. OOE+00 0. 00E+00 
6 O. QOE+00 QO, Q0E+00 








v4 0. B9E+04 
8 0. 10E+05 
9 0. 235E+04 
10 Q,Q0E+00 
it QO. Q0E+00 
ee O.10£+035 
13 0. 10E+05 
14 0. 235E+04 
15 QO. 90E+00 
16 0, OOE+00 
1/7 QO. 7BE+04 
18 QO.10E+05 
ey QO. 7B8E-O3 
as QO, O0E+00 
“ QO. QOE+00 
PiaP Q.Q0E+00 
Pua QO, QOE+00 
Pas QO. Q0E+00 
Zu QO. 00E+00 


THE GRID LOAD FACTOR = 0.186702 


MAXIMUM MOMENTS AND LOCATION ON 


RETWEEN NODES MAGNITUDE 
HORIZONTAL MEMBERS 
7 AND 8 O.1170E+05 
<« AND 12 QO.1222E+05 
VERTICAL MEMRERS 
7 AND 12 0.1222E+05 
8B AND 13 QO, 1222E+05 


** INDICATES + YIELD MOMENT 


*%e#*## ITERATION NUMBER [5S a RRKHE 


ADDITIONAL CONSTRAINTS USED 


QO. 190£+05 
QO. 10E+05 
O. 6/7E+04 
QO. O0E+00 
O. Q0E+00 
QO. 10E+05 
0.10E+05 
Q, 10E+0OS 
0, Q0E+00 
Q. O00E+00 
QO. 22E+04 
-.11&+04 
QO. 10E+05 
Q 


LOADED MEMBERS 


LOCATION 


32 620 


2. 000 


7 
~—s 


% HE 
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AUTOMATIC CONSTRAINTS WERE USED 


RETWEEN NODES LOCATION TYPE USED 
7 AND 12 7-O INEQUALITY 
B AND 123 3. O INEQUALITY 
12 AND 12 2.0 INEQUALITY 
MOMENTS AND LOAD FACTOR 
NODE NUMBER HORIZ MNT VERT MNT 
1 0. O00E+00 0. 90E+00 
Pa GO. O0E+00 0, Q0E+00 
i 0, QOE+00 QO. Q0E+00 
4 0.90E+00 0. QO0E+00 
ou 0.900E+900 0. 90E+00 
& 0, O0E+00 0.900E+00 
eA 0. 99E+O4 0. 62E+04 
3 0. 56E+04 0.89E+04 
9 —.17E+04 0. 10E+05 
10 0. 900E+00 0.990E+00 
11 QO. 90E+00 O.O90E+00 
Ze 0.61E+04 0. 99E+O04 
13 0. 19E+05 O. 71E+04 
14 0. 90E+04 QO. 10E+05 
15 0.00E+00 0. O00E+00 
16 0. O00E+00 0O.00E+00 
Pi 0. 1OE+05 ~. 68E+02 
18 0. 10E+05 —. 14£+04 
19 meee ale 0. 10E+05 
20 0.900E+00 0. 90E+090 
a 0,90E+900 0, QOE+00 
eee 0. 90E+00 0. 90E+00 
vias 0. 00E+00 0. 00E+00 
24 0. 90E+00 0. 90E+90 
20 0.90E+00 0. Q90E+90 


THE GRID LOAD FACTOR = 9.16E+0S8 


MAXIMUM MOMENTS AND LOCATION ON LOADED MEMBERS 


RETWEEN NODES MAGNITUDE LOCATION 


HORIZONTAL MEMBERS 








7 AND 8 

12 AND 12 
VERTICAL MEMBERS 

7 AND 12 

B AND 12 


** INDICATES > YIELD 
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QO. 1OS58E+05 
QO. 1000E+05 


OO 


© LOQOQE+O5 
- LOOQE+OS 


MOMENT 


Fee 2 | 8 
3. 080 


3- O00 
3 =. @@0 
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EXAMPLE 3 


GRID CHARACTERISTICS AND LOADING 


NUMBER OF NODES HORIZONTAL = rm) HORIZONTAL LENGTH = 190.00 


NUMBER OF NODES VERTICAL = =) VERTICAL LENGTH 


ROUNDARY CONDITIONS 


TOP CLAMPED 
BOTTOM CLAMPED 
LEFT CLAMPED 
RIGHT CLAMPED 


HORZ BEAM PLASTIC RENDING MOMENT = ©O.190E+05 
VERT REAM PLASTIC BENDING MOMENT = QO, 19E+05 
ZXSLP ERROR NUMBER = QO 
POINT LOADS ADDED TO THE GRID 
REF NODE DIR FM NODE DISTANCE FROM NODE MAGNITUDE 


7 ON THE NODE 0.900 1.00 


LINE LOADS 


BETWEEN NODES REF NODE MAG OTHER NODE MAG 


190.900 


LOAD NUMBER 
1 


LINE LOAD NUMBER 





pea) ta 18 


NODE LOADS USED IN CALCULATIONS 


NODE NUMBER NODE LOAD 


1 Q.9Q 
Pa 0.Q0Q 
mS QO.00 
4 Q, 00 
o 0.00 
6 0.Q0Q 
V4 1.00 
3 0.00 
9 0.90 
10 0O.Q0 
11 0,00 
Fa 0.00 
1.2 0O.Q0 
14 
16 0.90 
17 Q,Q0 
18 0.00 
19 0. Q0 
2Q 0,00 
fe Q,0Q 
ee Q,Q00 
oo 0,00 
24 Q,Q0 
2 0.00 


eeee% TTERATION NUMBER IS 1 HHH 


ADDITIONAL CONSTRAINTS USED 


NO ADDITIONAL CONSTRAINTS WERE USED 


RETWEEN NODES LOCATION TYPE USED 

MOMENTS AND LOAD FACTOR 

NODE NUMBER HORIZ MNT VERT MNT 
1 QO. QOE+00 QO. QQE+00 
2 OQ. QQE+00 —-. 1QE+05 
= QO, QOOE+00 —. LOE+OS 
+ 0, QOOE+OQ0 —-,. LOE+OS 
oO QO. Q0E+00 OQ, 00E+00 
& -. 1OE+05 QO. QO0E+O00 
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O.190E+05 O.10E+05 
—. L0E+05 0. 80E+04 
—. LOE+0OS —-. 10E+95 
—. LBEFOS 0. QOE+00 
—. LOE+05 Q.00E+00 
O., 20E+04 —. LOE ROS 
rs =. LOESOS 0. 50E+04 
14 —. 80E+04 —~. 1OE+905 
15 —.1O0E+0S5 QO. Q0E+00 
16 —-. 1OE+05 0. QO0E+0N0 
17 =. 1GEf05 —. S50OE+04 
18 —. 20E+04 —-. 1LOE+05 
Mw —, ,0R+GS —. 50E+04 
Pe —. LOE+05 0. O0E+00 
21 0. O0E+00 0.Q00E+00 
ae QO. Q0E+900 —. 1OE+05 
oe QO. Q0E+00 —-. 1OE+05 
= QO. O0E+00 =, KOE+OS 
Pega) 0.90E+00 0. QO0E+00 


THE GRID LOAD FACTOR = 0.80E+04 


MAXIMUM MOMENTS AND LOCATION ON LOADED MEMBERS 


BETWEEN NODES MAGNITUDE LQCATION 
HORIZONTAL MEMBERS 


VERTICAL MEMRERS 


*#* INDICATES > YIELD MOMENT 
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EXAMPLE 4 


GRID CHARACTERISTICS AND LOADING 


NUMBER OF NODES HORIZONTAL = 4 HORIZONTAL LENGTH = 10.00 
NUMBER OF NODES VERTICAL = 4 VERTICAL LENGTH = 1090.00 
ROUNDARY CONDITIONS 
TOP SIMPLE SUPORTED 
ROTTOM SIMPLE SUPORTED 
wer | SIMPLE SUPORTED 
RIGHT SIMPLE SUPORTED 
HORZ BEAM PLASTIC BENDING MOMENT = 0. 1O0E+05 
VERT REAM PLASTIC BENDING MOMENT = QO, 10E+0S5 
ZX=LP ERROR NUMBER = Q 
POINT LOADS ADDED TO THE GRID 
REF NODE DIR FM NODE DISTANCE FROM NODE MAGNITUDE LOAD NUMBER 
LINE LOADS 
RETWEEN NODES REF NODE MAG QTHER NODE MAG LINE LOAD NUMBER 
= AND 6 1.00 1.00 1 
6&6 AND 10 1.00 1.00 Zz 
2 AND 6 1.00 1.Q0 2 
& AND 7 1.00 1.00 t 
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NODE LOADS USED IN CALCULATIONS 


NODE NUMBER NODE LOAD 


x#2*2% ITERATION NUMBER IS 


1 0.Q0 
Pes 3200 
> 0.00 
+ 0,00 
J 7.00 
& 29.20 
7 J. O00 
8 0.00 
7 0.00 
10 2 OOD 
11 0.00 
= Q.00 
ms 0,00 
14 0,00 
15 0.00 
15 QO, 00 


1 HEHEHE 


ADDITIONAL CONSTRAINTS USED 


AUTOMATIC CONSTRAINTS WERE USED 


BETWEEN NODES 


LOCATION 


MOMENTS AND LOAD FACTOR 


NODE NUMBER 


pat 


QOD NO OD A he 


HORIZ MNT 


QO. O0E+00 
0. OO0E+00 
QO.Q0E+00 
0. O0E+00 
0. QQ0E+00 
O. 10E+05 
ge 2) ew Os 
0. Q0E+00 
O. QOE+00 
0. 1Q08+05 


TYPE USED 


VERT MNT 


QO. QNOE+00 
0. OOE+00 
Q, QOE+00 
0. Q0E+00 
Q. QOE+00 
QO. 1O0&8&+05 
O.L1OE+05 
0. O0O0E+00 
QO, Q0E+00 
QO. 45E+04 





11 


e, 
—_ 


1z 
14 
15 
16 
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QO. 190E+05 
O, OOE +00 
O. Q0E+00 
O. OOE+00 
QO. Q0E+00 
QO. QOE+00 


= 78BE-OS 
0. QOE+00 
Q. O00E+00 
0, Q0E+00 
QO. GOE+00 
0. Q0E+00 


THE GRID LOAD FACTOR = O.18E+03 


MAXIMUM MOMENTS AND LOCATION ON LOADED MEMBERS 


BETWEEN NODES MAGNITUDE LOCATION 
HORIZONTAL MEMBERS 

2 AND 6 QO. 1QQ0E+035 19.00 

& AND 7 9998. QO. LQOQOQE-O1 
VERTICAL MEMBERS 

<2 AND 6 O,. LOQOE+OS 19,900 

6& AND 10 0. 1O34E+05 1.740 HH 


**% INDICATES > YIELD MOMENT 


#*** TTERATION NUMBER [5 LL MEKKE 


ADDITIONAL CONSTRAINTS USED 


AUTOMATIC CONSTRAINTS WERE USED 
BETWEEN NODES LOCATION TYPE USED 
~ AND 6 10. INEQUALITY 
6 AND 10 1.9 INEQUALITY 


MOMENTS AND LOAD FACTOR 


NODE NUMBER HORIZ MNT VERT MNT 
1 QO. QOOE+00 QO. QO0E+00 
Fis QO. QO0E+00 QO. OO0E+00 





3 QO. QOE+00 
4 QO. Q0E+00 
ra Q.QQE+00 
& QO. 10E+05 
7 QO.126E+04 
8 QO, Q0E+00 
9 QO. QOE+00 
10 0, 10E&+05 
11 0. 65E+04 
12 QO, OOE+00 
15 0. Q0E+00 
14 QO. GQOE+00 
15 0. 00E+00 
16 0. QOE+00 


THE GRID LOAD FACTOR = 0.18E+03 


MAXIMUM MOMENTS AND LOCATION ON LOADED MEMRERS 


RETWEEN NODES MAGNITUDE 
HORIZONTAL MEMBERS 
2 AND 6 QO. 1QO00E+05 
& AND 7 QO. 1QQ0E+05 
VERTICAL MEMBERS 
< AND 4 0, 1OQQE+OS 
6& AND 10 QO. 1O00E+05 


*#* INDICATES > YIELD MOMENT 


0. 27E+04 
QO. 35E+04 
Q,. QOE+00 
0,.Q0E+00 
Q. QOE+00 
QO. O0E+00 
QO, ODOE+00 


LOCATION 








es 


EXAMPLE 5 


GRID CHARACTERISTICS AND LOADING 


NUMBER OF NODES HORIZONTAL = 5 HORIZONTAL LENGTH = 190.900 
NUMBER OF NODES VERTICAL = 35 VERTICAL LENGTH = 19.00 


BOUNDARY CONDITIONS 


TOP SIMPLE SUPORTED 
ROTTOM SIMPLE SUPORTED 
LEFT SIMPLE SUPORTED 
RIGHT SIMPLE SUPORTED 


~O. 


HORZ BEAM PLASTIC BENDING MOMENT 


20. 


VERT BEAM PLASTIC BENDING MOMENT 


ZXSLP ERROR NUMBER = @) 


POINT LOADS ADDED TO THE GRID 


REF NODE DIR FM NODE DISTANCE FROM NODE MAGNITUDE LOAD NUMBER 


LINE LOADS 


BETWEEN NODES REF NODE MAG OTHER NODE MAG LINE LOAD NUMBER 


4 AND re 1.00 1.Q0 1 





=a 


NODE LOADS USED IN CALCULATIONS 


NODE NUMBER NODE LOAD 


##*#% ITERATION NUMBER IS 


ADDITIONAL CONSTRAINTS 


AUTOMATIC CONSTRAINTS & 


BETWEEN NODES LOC 


MOMENTS AND LOAD F 


NODE NUMBER H 


DONG UAHNM EH 


THE GRID LOAD FACTO 


MAXIMUM MOMENTS AND LOCATION ON LOADED MEMBERS 


1 0.00 
Fa 0.00 
5 0.00 
& 3200 
a ver os 1 8) 
& 0.90 
7 0.00 
8 0.00 
9? 0.90 


1 *#*## x 
USED 
FRE USED 
ATION TYPE USED 
ACTOR 
ORIZ MNT VERT MNT 
0, QOE+00 QO. O0E+090 
O, OOE+00 0. QOOE+00 
0, 90E+00 0. Q0E+00 
O.O0E+00 O, QO0OE+00 
QO. 20E+02 « -OE+02 
0, O0E+00 0. Q90E+00 
Q. QOE+00 0, Q0E+00 
0. OOE+00 O. OOE+00 
0, QOOE+00 0. OOE+00 
R= 2.4 








meme Se io oe 


RETWEEN NODES MAGNITUDE LOCATION 
HORIZONTAL MEMBERS 

4 AND a 46.87 6.290 * 
VERTICAL MEMBERS 


**% INDICATES > YIELD MOMENT 


#**#% ITERATION NUMBER IS a RRHKE 


ADDITIONAL CONSTRAINTS USED 
AUTOMATIC CONSTRAINTS WERE USED 


RETWEEN NODES LOCATION TYPE USED 
4 AND ra 6.5 INEQUALITY 


MOMENTS AND LOAD FACTOR 


NODE NUMBER HORIZ MNT VERT MNT 
1 0. O00E+00 0, OOEF+00 
Fis Q. O0E+00 QO, QO0E+00 
ms 0, OOE+00 0, Q0E+00 
ae 0. OOF +00 Q.O00E+00 
oa 0.156+02 9. S0E+02 
& OQ. O0E+00 0. OOE+00 
v4 0, O0OE+00 OQ. 90E+00 
8 0. O0E+00 QO. QOE+00 
? 0. OOE+00 0. O0E+00 


THE GRID LOAD FACTOR = 1.8 


MAXIMUM MOMENTS AND LOCATION ON LOADED MEMBERS 


RETWEEN NODES MAGNITUDE LOCATION 
HORIZONTAL MEMBERS 
4 AND oa £0.00 6.250 


VERTICAL MEMBERS 


%* INDICATES > YIELD MOMENT 








i 












































































































































































































































































































































































































































































bw a a oa.” & ' ——_ ." 4 . ee oe cl er 
f Sf “h i $ . o-3 8 ot Rig 4 t 
4 4 a, : 
Vit O'R b t ‘ oso 4 ot 
s " at as 
TM MAG vs Saas Pe % t - 
et aes War 
eaten + iq 4 we thesH296283 r an) 4 4 ery 
Po SENS teres Gti he ‘a ; ' ' > » ' 
Ce is m Limit anal ft ly loaded eae 
Ropeeemeey § Limit analysis of transversely loaded gr J | | 
ted : wy f ‘ : Fi ' 4 
: : 4 1 008. ] WI 1] { ' y Rn Th) { 7 u 4 ) 1 
‘ * kd b Ae tk p I ] il] | ' Pir ' rs 
. ome et, : , ey ; / J * é . 2 
wor ad baie ee ha € ; . hea tira os ot i { | | re wae : » § ‘ o 
PAP ER: ony Leia ody : ah Nes hy by j , | | } | ° : to . . 
: Pt ite FE %. 4 = i as 5 w6 0 ch ay . H Hee | | } j ae 1 4 , ’ 
4 Ea hye att , Wh | | ene rT n : ; hie : : 
Aut. A hee ns KE! | a ‘ P a! ‘> : F 
: ry Sa nie a &% +4 aoe : 4 af a | | | | iT | 1 1 | a e ‘ 55 or : 
3 OF hey Sage robe 3 Aare" ed 4 _ ' 
oe tae wb a hee ey Pokey fia eerie ‘ ‘ ' : 
ACD YUAN KTP 3 2768 000 64769 7 oe lee eee ee aL 
bo OA 8 ol, Rt ett, FObSib ..* suNeae © 4 4 soe hee) ae : 
yes rs er ier au é be s : $ ® 
eRe ae ea eae DUDLEY KNOX LIBRARY ‘ey pee ; , | 
SG Enh ae, ry miyh 8 Wig cw } ms {4m ° rs 4 
M bdnlwe pose Se ‘ b Unda "a4 Pace to us : "fh ofa 
Ooh ots edt ak op . ! ‘ 1 uJ La: a "ag ores aos ‘ t . Fe - a - ' a ‘ - s f ' 5 
| et ie ae ~ , t Ade 4 roy L F ; . vol an , 5 4 : ; 
PNT RE ot ot wt 04 ot, ig , S aan rs LA tok © eee oe ee ah A a a! * ; ' ¢ Li , ' ‘ : m4 5 
1 OAR oe tery Sapo : , ‘ tan tO “ ‘ ; j j ne 4 ‘ ‘ ‘ : 144 1 sa ' 
*} re ter 1 , rege & + bene t Pi 4 ¢ 4 c 1 ® ' ‘ 
va 6 O14 de #176 1B Yh OLA AH OU Cae Oran re > H : » Cer ceed i, i | “phar Ve ¢ D ' 
Pita bbl AL St Ar of ol het wh tg mee g gy 6. toe ' Regu fon : A ' 
re 0 IEA 2 8 eNews PAL 6 OO LS a 4" o3'4 Aa ‘ LA .! oe LO of ta ae y a ' 1 P . a t 4 
ret Sle tsa hb. Ite eng? Od tye Mel 4 ¢ yan : 3 1 i Pere A ‘ ; . : vs ; 
i el Cet a eS eraean Pe OS ae es Oh igs SLY Car) a’, a ¢ re i te Sa ' 
Wirt Pw rhed) Baan Caen tit Oe tls fee Sat a Vn gy OUST LISSA Tig ies us ' Inare : L 
adm FA SH Keer oy Banal the ne ig Oa F a Lay Ve ofa ‘ t ‘ ' - ‘ ' 1 ’ ' 1 1 
J divans Oey erers yr we prt Sieg A Mert i ghee wh Wome 0 tts th ie ah aad ae ‘ : ; 
> "h. 2 . - : boat 1 & 1 ' 
SHA A ht oh i ‘ ‘aa if th 0 Ci ait ast 4e 4 Pay : 14 ° ; 
se dk niabie de ie El bie ly 4 dt ct He ; Meese : 3 mt ‘ tf ) : ; Ae : po 8 a " elg x ; ' ! re 
mi e . 4 « 4 a ¢ 4 J 
Yee HE ane AHA Sm Bm 4s , 4 me " ‘ . Ad 4 ee . 7 ' : ae | ‘ys = , 
‘ MAE one Pe ore a z x Veqe yt ok? epee ' ; : 
i De By Das . 4 ¢@ Ln . 1 
A Mats a ME ay “ ee warren wis ata eu ; ie ! 2 Se a age X ; ; , : 
| alae cate Cob Ll Oe 0 ond S.ep ne weed ci reer ER US? Ce Bras &2%41@ 8 : a} are A Ln . F ' : ® 2 ae 
Be oti FAOMEA A DS be 205 uns apy Se faye ane ee aM owen aid steiGar ae enna oF i aus j , ‘1 : 
, H Bit I Md ee ; ek Oe wi. f fuer 1 i 4 by . 
ad! Cte at SfItnm Aa Leenga A ! (een 1 a 
Sits hve s rece e 4 ACT ICTS (Mi A eade ee Ja bie Sas ! ' ‘fC wn se ie j 
SUN A Meat Pas Aceves MU Otret eas’ he ee a ace eae | ' , ' 
roe Oty yh fin eye EPR id 88 hcg roared rie Pe er ae ak Oe ater nies i 1 - : wee 
We : : age ‘ a : a¢ 
Pe aie) late a let OOF F616 8 wth og ah . A ius ads ea" 8 . : ' 1 F as =f ALIN on t “1 vas te Hq ' 
Dr eh ee te lak Penne en Ae a ae fides Poteet OIG oe tM, fee ‘ aicT eh tage ; praeeee 4 t 
o's 7 sh as » sek * 1 if ' 
3 OP eb hid Go yi gate nedweduws oy as tee Mess eee sy pots 21) a 5 ma 4 ‘ Paw sel é pr <"" ‘ ' E e ae a : er ‘ 1 rs Py ; 
Saw OU) eh Bhs e on a eM oe re ais 6 a AS FAH RA OO ate Ot OA SMH IIe 2 re taf 44 er uye ets esses ' : oo. ' 
ae agers OOP Sieg teewe a cited be de egl@ is ares Stott ebal “4 aot nes eee Lt e os FO ues j ‘ ao atrayn “oe : % i 5 ‘ . 1 . ec Ch ms : ta <i a Le 
FOR etl Fit oh ee we Pres dat Fs R oer 2 8 .bas , ' Ce ey os re ia 
é rs Wate a oe so18 oe aS - Gass iS # « « 1 
tan aoe re te I ON minis geared ME hog ed Se Fe ae: f2 treme Ur ee Ae lee Amr cee aus fy at ts 3 ne : 
Poe tag 00, sa CP ek Per ee ale OP E04 bw Heh ete enw es a J ay ay, ' ; ' a ie seul te as fe a ; ; 
0 Pe aes. et EM abil NW PL eg lt! to bt bed 8k 04 3 . 7 P PT CN 1 ' Ges ' aos 2 4 ‘ : 
em MF ad iat Moor ts pF PO pod 5 '‘ 0 bh beaet “ans, wae fh oem : q' 1 
Ne EN formed seo a Oye Out ig ab Caigs s That se haa 4 “6 paenge cuem as 1 fw éotseaneas oa ‘ . 4 : ' ra 
Pie Re cau WOH ne One Ol se rte dat of 40k hte cae * oF ePrn 6 0 of d ode 1 tu ‘a . ‘ ' 
nh Asi AMee LD eet eas MP ihe ase tg tA Oh deg ; a : 
ee ot Contd 6q- tip iets Pet bar tt win RetegeL a UNC 27 ar ae Og Oe gag r Po 1 
ain gd nee Sohne Hise” 8, EEE 8 Pe Bn Caer oe ifade ‘ t 1 : te 
one © oe w 2 wet e 5 2 ot ob on i,’ ein ‘ : [nl 
wf ote ol, ey % - “Aw ons od @ue 1 ‘ ' ' ' ry 4 
SW sg etree ge 44 of Set A eaten Oe et eed “fied ati I : : : ' - vyer . “ 
PP PR | Ai OI ce Peg ieee Oi! Pl Ok gue kg tO OI kg 8 Ase : a? toe i ' 
PR eer eer e ate aawe. Be! ed ee tT io om o? 9 © ' erie A : a F) , , e 
fen . hd 1 SOFT Coe Be Oo . C4 4am Oa At e 4 
ahh, ee baud dk ie 3 cay Fay Oe oF ofa os 0 Fy rad oe Lt Pari ee Pe "od om y ‘ ; ' : s 
ol oF eid Bai orl aalnd Fal Hs eB iis SUP Melee nag gut Ost et aye po Of A eh g gt ven rice fea Rs : : 1 = Aas . 
% ae a od ‘ee ee etd . ' @niat i "ae P 1% ‘4 « 
mA MH OK mag © ey aT te aa Sega ! ee ites Dias Se , é a mo ; ; " 
cy I; FOP PEGE or thes hee se OE AL ans op oft WB rh Ot eg es b # Pie yg 1 ; s ‘ ‘ 1 « 
Pets Geen Fo hh Po eieies Ok de Seee 14040 Me ans sf; oe ee) aKO ate aoe , 
Rta gi pence Ce are oe gents ad ae ee ined ee A Oye ei C4 Viigo 4 ks i : 
q tS Boat sae A wea, hana Pre ha ty om” persed dig te ie teaewaeg « va ns, 8 (avrag 4 4 . < te : . 
Sb bel 3, OeAoke toe WO We le 8 hid oP Dod Gica Peet ome go yy eit ; ‘ ; 
A hw Ae 4g Seuniss ong we OF Tea Ng Bt oO 4, reid Fe oan da ih ete Peete. | ee 4 é « 
nal AA ee bas 9 he gf NS oie Pst F< did SLOT Bee nethe a , , | 
FORM ews hse Aus ee TM EO ret erg Peat FRSt ea Se eM sete eC “3 P woke 
“0 FS os Ce dies ae rt ae Samer’gi ‘as om . b . ; 
me . * mies ‘ 2 Ou = ~ #9 aot 
Beh geg aS fF fet dems ays oe gata eget te Nisiars) p sicg poke se the , , ‘ : 
' s . 
' a ' « : 
re 
‘ ! 
5 of « ' 
' 
a slamwe 2 ay fet * at gt ers Ate £4 : ‘ 1 ' 
ale : < er ge auto s La 5 fs s ' 
ihe cae EOS Bees rN Y Fs Mie Crs rating Aes Sec6@ Ma Me vestgr ot, ce a «. F F ‘ ' e « - 
at oF antic a caueee” Fah Fah 3 Py aide th f On *y wets Me Pye mane pret ne 
” tO fee sins bi ste Pudiwys 0 33 Apa ele: 4eyag - aie ‘ . 5 
Pers TV se oo oy dae * Peat Mi ‘raie's a 
tec OY Sad beragiope fet BEEN tage +t 
A Set cil a ILE a aed awa 1 ; 
a« ae | va ¢ es F) » yee e e bs 1 ' 
' 
. os Le) ee . Se 
Sep cwsuw we, - oe “au id Pi 
att Fer het, 89 RE tte : 
oat > fury 0 oro a é 
eye, ge ‘ 
' F . ' r 
: “ 
‘ . * 
' ' re ' 
' bd « 
° ' 
One = « 
* Py “ ‘ . ' 
‘ ' 
ioe . ; « aD 
1 og ’ 
. ‘ 4 i ; 
' : 1 1 
1 rs os 1 
' : « F = a 
¢ ' J . 4 A 
‘¢ . 
‘ : ¥ 
. 
7 ' 
: 1 ' «* : 4 af 
, : . N M : 
' ace ae ' 
. ! ‘ ' 
« 1 e 
Ld] « 
ie 
- a é F 
« ' 
' 
' * 2 5 : « 
. @. 
. ' 
~~ 7 >’ : er Uaie 
19° . 5 « 4 5 | 
OU he b . 
vn . ry 
. . 
vm ty % . u a a ° 
. . 
Sige «4 “os : eee 
f i ' 
awe 4 - ee er ; : ° 
ee a in er ' 
ph RO x ce 
a Sor ‘ Wate ty, 9, et ' 5 ' ; ‘ 
ar Ne Ma shes . i 
= > ee AeA . ar g i : ' 
~~ t < 
d ‘ - = May tas — a 1 
i Se Weems. wy heen es ee wR GED - Ee. . oy 
.. me ae PETAR, ated! AY > aw 
at 14> ogi me ay viens bias ex Peay CE dy . “VGA NOK, ms a F = é 1 
oh & 4 on OC d* RBs Meh y ives . WD? ostsm rm ~~ we «2% P =. FL tas aFx . ‘ 8 « ' 
L ete 2 F Y FF 38%, 6 : ~ - ayy >a i x | ' . ’ . ‘ . « 
: a arate a TINE NS ey *. mS CS ee LEC tate a ate ae ot ' oe ° ig ae ; ; 
z a . “tou Se br ate 3, « AW cere De ze is ‘ a ; "eae, va : . . Pt * H ’ « & q ® e e ° ar s, @ . 
) Fe ns r A \ i<_3 Mia fae ie ¥ e “Ww ve = ° 
Sy eet Ue VO acs Git ee) % Svyteres re Bafa Gee AEN Ae ay UN atiaaehs =” : . ge . % Ake - . ° : Pe « ’ on 
- - . aa ® ws Mw. My Gein ee " < ba % 2 Win « me . a . - t 2 7% oe 
% by 2 Vs en wi we ONE Qa ae tay ¥ t : « ; me Ca ’ , ' i e ‘ 1 
© Fr GA dem ma Shab en ae a eine ere wh hbk ay Face, 7 eo Oa : : ‘ : a : 
; WA. Se ron dee % MOM in At hy ELAS ht Sete way a. we Say 7 . ' ® 
2 We pR eek te ee! PSone ye NT Ne wate ETE EE, Bly e- soa ‘ ot , ‘ 
Nes Whee es te ET ie totale he Oe ee Gate ke ‘ ‘ . ae on ‘ 1 i . 
ba J Nee 1m, Ma SESS. a tdorary cnaie OR AN ey NT ee TENGE OO : oe : wt t as es 2 
eS WRG BN Ee fa ae PENI ae tee he Ore ree SN EUR a, Oh es as ~ oa « ° a ose ; 
= SUED ew + Noe eBere . ONT mar Moe habda a Pt) ea e ' a o« 2 2 ee 
- $N ~ FOFEF. arias . WS ar " mee Fi he 
Wee Ras, oy S% Ma Teme GY ee eins eee 46g an 3 PR DR ve a e « aie = s 
he RE Ye, ~ 2% wi tery aoe = pote Se | se Mp lon pa ee Serres he % 8 : P ; c : s 
Tees ha TE Re Kg ahre eee ech sat bade . ee FRR e KAM tig ee . . se . « ry * 
at eee tas wines COMA SR in * = pre ve a Sena Mala him te RTM otk EY wire . a ¥ a * « i 
aes x print et eee poe ay ee eo ay SYA Re we Ne ae, a TH ara be tte : \ a : ‘ ~ * & , 
be a tk me ES FUG - Fo my “ eet eee Ga ® ie Te eee a ‘ ot 5 " " 
mA Sm tarig amet, BNO FN to eT ONS he ee EE AN os WN : = sow Wey ' . 
. 1 Fa Nee Mog pag SU Ue ep TA ren sie Stoo we OIE Bw A Pr ng ne A oe . 2 gee ‘ ' Py . oy 
~~ "yh : = f x v Reoew B bs = 1: Li y . 
oat TAI Rerde t =e Ue RN FUR Lysine 00 ek Vers “OVA, U - .. z + s *. % 
SNS ee magna Ne “Lame Pata ney, wa We eg Rah ie 72% WI i « : ' ' F . 
Fie oder, CURRAN Mae Gin el ty 8 ce dy pict Ulead ee ee wre MAL Ber ary ¢ « 1 het | . ' aoe bs . 
@ ese Ne aR, Wen Fae “ns, N89 ; ; blah h a a he : a ' ' " 
ph ae halal | Th ty > % aan yey 3 ‘om pie: a, +" oe “ o STUN e a g, : Sh « a = s¢ ‘ e re i : f e fe tue 
‘ Sr tareele tote GT sd de P na sn WP Sangred oer MNS Cd a Paria ctr S : De Pa ede goed ene — ee 8 a : aete i omerars ‘ 
Ren a eee i nes At ip Colt 1 eh A Te oa ee ss ¥ “bo eo tie ff oie . ? wo ihe . . ' 
, Neate +a wER Nn eal ta BURUWEN G wit tw Teg ae et aus redial a Se Yt ¢ 5 P : « ae a Pn 
pan = Rage jana a ov. ms “a "Sig < fq . ae PeUt eA er ‘e a 1 1 A 
> mm caplet la ah gy a ee te 8 Pept Mop teat lk Reruns R Easy ween Be S elg  Aeeae ities ; a * et y . : ; 
SE eae Se fanin ad te tag eee RAR CY owe wis s pita) re hes Sh ee vee: vO % wy fee yeaa ° Pa yr eta le af we F 1 ‘ 4 7 =e 1 « * i i ' a 
"Woe "Eas eta pine tt he het * Pett we \eae Ven mer vn hilt Mee sae Pe ve te eve aw Oy tee’ Sane , erate i. » 4 : =e 4 
trees ees AO My ial teen “ EAA oe DH WG teee-y Fe NU Ok Serpe tiene Cute , é SCT Ta Ci : aes + ‘ ' Phe q : . ‘ Sr ataye ' ' 
tet ke tet MMM 1SUl ery ay @ igi ede LY Umy Qig Ke Aah 0B Weare moe ue git aftr Sf ee aa met aly y 1 es "| 3 r) F : 
PEs ee te : een gare x see . “ rn. te © a é er 1" 
“PPS ay CES Wi tel aia tals be Be SL ¥ pala p O'¥ethy tee ce bobo ‘e eee “hed Teavy ca; AS ORO a rice 1 ‘ Ws ; by P ae : sn ' ' ‘ 
Sa EN ru where eee y ; ‘ rk eee Lt ew et ay Oe. ioe: « et ate : ae = 1 Atty . epee 7 2 1 : ce ei ‘ . 
2 NWS Melt u's. 8) 915k 4 ve ant . , . VU UTA Rey mee _ 6 % ; aoe aes : = Cie ag . 
SRC ed tbe er pen AVANTE "we a tlh y Leah a yA sa R HK hal Le ne 4a é x : 7 ' : =. x ; as « « . ; ‘y . 
ne ay a he TT ET at imc ‘ Sara ae Sivas % Ne th, lle enn ha Cea ta Te Nl neater lee Ge rere aoe ‘ wa ce Sg: ; oe : 
A ww EUW, He PRUs Wana “O) SU my 5 Pay iver & PEO Ney en ty We i: se . vy tA ‘ eis : ¥ * * 4 
et : . we = ; eC wre i yea ’ yas H&E Y ] te wed r i . f 
NS NES BOA See ISNT CNwoe ee ee CT, ReV use ates wD VAY Gt 8 meg bt pO iT ie EU ih elieal 5 Cy aan on ts aie ence he wa of ° ; 
Ma are GN ae eR eR cm de ec ee eT AIR ona ae ee UES eth te ee ee UNO WOH De keh de he eee See ’ ; wen ahora tu 8 \ aus : ‘ ‘ 
ee ee ree er nt WO fu Weve ony | Steele ta BE MV ae Gh s,s We wee kN ues 1 eee ne ’ oe ‘ ; es . ¢ - ' : « 
8 AE ew ey A toe ee eee Be tk le Te Re OEY oY RR ere gt A ore te e ‘ e1 ‘ 
ee mace erate CE UN LAENY Ua Qiheecee Wye hs aan ahi MeN A wun a A we ee Ss: cuceryecteee ime eine : : ’ . a! wey 
5 ae ee Lod on. 7 ; ‘ See ; Fle Ge yD, Wee 8 ye we be, | ms wer ae e v wha . A = 1 a ¥ : 
BNE Sorte ay wt Bre &s ‘ee HE sb VeeS PRS Ke wg : re, « ee de ae ie " * t . ” v- a % o8 : + tas 
S/R mae ees We og, ‘ae SREY WLR EIS eae eA oH, owl Gy Va ay eergt ya u at ad ws » we LE a 1 . me a ; , { < ‘ 7 
3 ewe vee 1, Awe ye + 3 a ee re $ ony ‘ 1 ee en 4 1 ' 1p 
LS brace a ag waa ea kare race hae area eee Bier X& &? re OTe & mh a rence ahi és cCVUVE Veep . FS oie A ; : a : a A a % a e&-« : se Sa : : = 
yes 2 BOE Ae rg Win tees FR Ass We ree y HAWS Wy egy oa © : ven : ie Pde Ot harry ’ UJ : i L ' ee ‘ So 
Ne Sr me en eT Ae be Mae Mbiee rtdede oa he ta CHAD Wray Pits de oe hh Se oe Pe, LWewt we wee ' pee os UL veer ay ea 4 . ‘ Leer et Jar 
wee Oy me Few Te Bee oy wy! p ie oe BUR Gwar SEEN & © beh Ue h ee Ae eee ee “4 ir | + 64 *o4 $ 4 « 6 ' « , 
hee Qe BH, gs on f ‘Sie-etWin BE Qe Se Syme, Whe Lom tas hina th aL ae ee ‘ ‘ wy we ‘ ‘ys as 5 
- te ot ef Wee Bite ne Wy payee HG Kg W Wray yt rd: St BEI ogy P a 27S OTW 4 2 . be 
VW RUE Sra, AU SV ene Oe E OE HEN SOM Ye ane: 8° CYR eo Airh td WOU 0-58 go VOtON% « (rpen @4s vi a , ¥ a 1 “os Hur ‘ « y . 
.) Awe a eee, SEA HS Ora ee ere Srv e'N Gey, Samu: aoe C2 weet & Lew ee dada te te hs eee Ee UU WS Ue ai eien Vane! ber La Ie : ' : -— ‘ ne . : 
We Ulery, “Sere Oy oie args Ure & hy 27 Vee Req & 0) ,t SS Eh te Wey S005 hte LETHE wee 4 awe be fof Lt ere : or Pacha La a % c/o e's eV Va. % ee 5 i : . “i . eu 
NI 9 tt ea haven eT Ee eh ich od = ry me yo etl ARE RY Uren pb adie TAT %e rere re tah HO ik ee +8 a ‘feet Wy at Sisiet es ee i , ‘ om 
1 trhealan ton! rae ne ae, a een , “ath? beth? za. ¥ Oe Hae Ee ene te i Bea Se OS oe Sear i Br aac fotiw ¢ +4 ae Eee . . . ; Soe e 
Weta bot. t. eee OG Fb ye WIG 6 Pure UTE Me Gat ee ET hay Oe b« : -§ owen : Ly Les SR t i ' ! : «4 ' 
NS, tee ewe FEY 4 SMe yO ue WW Oy 1S Wee oe wv PM ieee c ney eed vs ¥ BEA Ts ‘ « at . ; et, F ag . ‘ aa P ; 
A ek serets SOG ! : * aw uNa sae TY hale tet Yor oe Mel "UP e ws re i ts! 4" trey ae! ; : 
" ew my tase rae ‘ : ews Ane ae, ie Pe Tee . y Lime ws r oe sa 4 he . cere ‘ ' 7 
0 Re sO tte ast aye hn EES Bok Vteces raed te TE WR CA Bred tyre ate et LOE NS RT 8 Ph Re , ‘kon apie ie roa “i i os Batre o. 8s : r 
vA Werk Fa ad edb temas ee k : Jee aia et ang 
Fase et " eons ve ‘ 1 mK a 
aR « a | - ¥ : Be digetige mite ‘ ber ' 1 . 
p &, ' ee eet ie | 
& late, We | . Wee 418 Ga 4 ot ered us at i) : a 1 ] "4 i e 
ieee AU ergean : * te : ; : ? 
jr ee ey neat Mreeusy Wega sane ve aun, haat = Nea : 1 : gona 2 R ues LS oa oe : ‘ ‘ a ‘ 
naka 7 . % 18 > 
. ar e wt e 
a § ~ 8 J t 7 - 
See ~~ 4 UNE) are ‘ eos : i a : fae 
Tt 1 . ’ 3 a - «* Ph e 1 « > 
SOL Tat ba ee eae ey EVN we Ore, ee : ee aes’, Pe run , Ss . : . : ' . eae 
BSAC OS La teehee ot Erte Chi SEL pin ORONO ST NTA : icay oe Knee eet e sais ce uers* : 
a Mb x é a a x "ewe ‘ ANY EY O62 be ’ : at ¥ J ‘ ‘ res “ ' 
Ae Fi hate eG Yawerie hi Sy 5 aye Re table fe 8 eREtna rere i | P : 
“eh Arment te ent etek & tim HU u Va We, ee pees tt Cl vee a el a o« g ' ; 
Wty OOO ea weg ee UREN. 8 Rh Gratin WEB HIE BEG A a et UU Le gine : vote 4 iets ' ‘ ; 
ee Ge RYE As Bree %, be 4 FE atte te 0 ae hwy . t 4 Soma ee ; ; 
4‘tncms otetat SAA MSs : e eee ure ee pe Lh meine SS es — 2 »4 .t wou 5 les mt : 
Pe eh Oh Ad ek RS ve 4 A Tree as ‘ : Stes ae 
> BAAS tena +e t vuminane : Mae gt = 1 on aes a td rae e 1 
b 3 t bo | ; 4 
ANS 19.8, 08'S Why, ‘sy he ae ‘4 $ . +o : 
ihren ee NS. Santer make h Wie tt in tow tites ; : ee . ; 
Ne, WM 4 oe BL em | th. bs Ue Wilt oven VE Cea Bete ‘ ¥ é Vou & .» Ls es «= Fea 7 A ' S ; 7 
ee RAS A OU te yas ce “WGIO ey Ub Oe A age Cs eee eee reac ‘ ” 
Ae Gd eh Owe ~ y & ; PingtesP se “Or.” ‘ 4 
ey ® . + €@.3, : TPR, May s,¢05 os rn 4 . z Pa ’ e ‘ ¢ 
BORAT CEU UTE 1 EO eae EI 4 mytete yee Ue OLY a wg Rea ar he ame taty Ue e *, ' ae ; 
heb Sis ta nth lt 1D a Aare vow hol 6 e am inn yy : CB Ss Tig oe eee ‘wan ale ae : 
y Da ie tee 4 e ts : i ’ ’ 
vie rere y. “ney it Wie ges ‘ wee pt! wv mews 190 | 1 3 ‘ " ' tals ' ‘ 
WY MWR Wene at fa ae he c oh Be a Wns \ ' . e o4 « 
VV. noite cays uta ' i? Or1ip gs . ’ ' ; 
Cw breee ee OEE t et Ne yetele TONED Db pg? OO 8 : ' ' 1 
NY HOWE cr he ni Sag 1°. bs as cul die ee ee eee = tr12#) . : « oe acs ‘ . 


